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We uncover a disorder-driven instability in the diffusive Fermi liquid phase of a class of many- 
fermion systems, indicative of a metal-insulator transition of first order type, which arises solely from 
the competition between quenched disorder and interparticle interactions. Our result is expected 
to be relevant for sufficiently strong disorder in d = 3 spatial dimensions. Specifically, we study 
a class of half-filled, Hubbard-like models for spinless fermions with (complex) random hopping 
and short-ranged interactions on bipartite lattices, in d > 2. In a given realization, the hopping 
disorder breaks time reversal invariance, but preserves the special "nesting" symmetry responsible 
for the charge density wave instability of the ballistic Fermi liquid. This disorder may arise, e.g., 
from the application of a random magnetic field to the otherwise clean model. We derive a low 
energy effective field theory description for this class of disordered, interacting fermion systems, 
which takes the form of a Finkel'stein non-linear sigma model (FNLcrM) [A. M. Finkel'stein, Zh. 
Eksp. Teor. Fiz. 84, 168 (1983), Sov. Phys. JETP 57, 97 (1983)]. We analyze the FNLcrM using 
a perturbative, one-loop renormalization group analysis controlled via an e-expansion in d = 2 + e 
dimensions. We find that, in d = 2 dimensions, the interactions destabilize the conducting phase 
known to exist in the disordered, non-interacting system. The metal-insulator transition that we 
identify in d > 2 dimensions (e > 0) occurs for disorder strengths of order e, and is therefore 
perturbatively accessible for t « 1. We emphasize that the disordered system has no localized 
phase in the absence of interactions, so that a localized phase, and the transition into it, can only 
appear due to the presence of the interactions. 

PACS numbers: 71.30.+h, 71.10.Fd, 72.15.Rn 



I. INTRODUCTION 

In many experimental situations, the effects of both 
static ('quenched') disorder and interparticle interactions 
may play comparatively important roles. These issues 
have, for example, come again to the forefront of de- 
bate in view of discussions centered around the fasci- 
nating, yet still controversial "metal-insulator transition" 
observed in 2D semiconductor inversion layers. 1 Unfortu- 
nately, theoretical descriptions of quantum many-particle 
systems incorporating both disorder and interactions are 
typically quite challenging, and difficult to analyze reli- 
ably. 

With this state of affairs in mind, we apply in 
this paper a powerful analytical technique, known 
as the Finkel'stein non-linear sigma model (FNLcrM) 
formalism^ to a model of spinless lattice fermions sub- 
ject simultaneously to both, static disorder and short- 
ranged interactions. The isolated effects of disorder or 
interactions upon the lattice model that we consider are 
already well-understood, yet still non-trivial. Our goal 
is to gain insight into the interplay possible between 
the Fermi liquid, Mott insulating, and what we term 
"Anderson-Mott" insulating phases of interacting many- 
body Fermi systems in spatial dimensions d > 2. The 
FNLcrM formalism admits a renormalization group anal- 
ysis controlled via an e-expansion in d = 2+e dimensions, 
permitting us to address the following general questions 
within the context of our model: in the simultaneous 



presence of both, disorder and interactions, (i) does a 
conducting phase occur in d = 2, and (ii) is there a 
metal-insulator transition (MIT) in dimensions d > 2? 
If yes, what is its nature? 

The primary purpose of this paper is to present a de- 
tailed derivation and a thorough discussion of results pre- 
viously announced briefly in Ref. 5. A short summary of 
these results and an outline of this work appears below 
in Sec. QZ3 



A. Interactions and sublattice symmetry 

We study a class of "Hubbard-like" models 6 for spin- 
less fermions at half filling on bipartite lattices, possess- 
ing short-ranged interparticle interactions and quenched 
disorder. We work in d > 2 spatial dimensions through- 
out. For concreteness, we consider in this introduction 
the hypercubic lattice. By definition, any bipartite lat- 
tice may be subdivided into two interpenetrating sublat- 
tices, which we will distinguish with the labels 'A' and 
'£>.' The d = 2 dimensional example of the square lattice 
is depicted in Fig. Q] 

Our starting point is the clean (zero-disorder), gener- 
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FIG. 1: The square lattice. Labels 'A' and 'B' denote sites 
belonging to the black and blank square sublattices, respec- 
tively. 



alized Hubbard-like Hamiltonian 

Hp = - t c Ai c Bj + H.c. + V £ Sfi A i5n B j 
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where and c Bj - are creation and annihilation opera- 
tors for spinless fermions on the A and B sublattices of 
the bipartite lattice, respectively. Here, i and j respec- 
tively index the A and B sublattice sites, and the sums 
on (ij) run over all nearest neighbor A-B lattice bonds, 
while the sums on ((ii')) and ((jj')) run over all next- 
nearest neighbor (same sublattice) pairs of sites. The 
homogeneous hopping amplitude t is taken to be real. 
The operators Sua/b = ( c a/b c a/b — §) d en °te devia- 
tions of the local sublattice (A or B) fermion densities 
from their value at half filling. Finally, the interaction 
strengths V and U appearing in Eq. (TT]) couple to near- 
est neighbor and next-nearest neighbor density-density 
interactions, respectively. 

The model at half-filling, given by Eq. (JTJ), possesses 
the following special symmetry, which we refer to here as 
"sublattice" symmetry (SLS) [this symmetry is termed 
"chiral" in the classification scheme of Ref. (see also 
Refs. i@IMillllllllllli3li)]: the Hubbard-like 
Hamiltonian in Eq. ([T]) is invariant under the transfor- 
mation 
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where we simultaneously complex conjugate all scalar 
terms in the Hamiltonian. [This transformation, like that 
of time-reversal, is antiunitary. In the presence of time- 
reversal invariance (TRI), SLS is equivalent to the usual 
particle-hole symmetry.] As is well known, the Fermi 
surface 19 of the half-filled, non-interacting model, Eq. (JTJ) 
with U = V = 0, possesses perfect "nesting," 
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where e(k) is the non-interacting bandstructure, and Kn 
is a nesting wavevector. For the hypercubic lattice with 
lattice spacing a = 1, Kn takes the form 
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where the numbers n s — ±1, with s E {1, ...,d}. 
Fig. depicts the Brillouin zone (BZ) associated with 
the square lattice shown in Fig. [1] The set of nesting 
wavevectors {Kn} defined in Eq. ^ span the sublattice 
Brillouin zone (sBZ), appropriate to the A and B sublat- 
tices of the composite bipartite lattice. For the special 
case of the square lattice, the boundary of the sBZ also 
serves as the Fermi line at half filling, shown in Fig. [51 

Nesting and SLS are tied together. Under the transfor- 
mation given by Eq. ([5]), the hopping part of the Hamil- 
tonian in Eq. (TT]) transforms as 

f ddk t f ddk * ( \ t 
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where Kn is a nesting wavevector as in Eq. (j4|), and 
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Fermi surface nesting is, in a sense, the defining 
property of Hubbard-like models for interacting lattice 
fermions in d > 2 dimensions. It is the nesting con- 
dition which makes the ballistic Fermi liquid phase at 
half filling in such models unstable to Mott insulating 
order in the presence of generic, arbitrarily weak in- 
terparticle interactions . 22 ' 23 ' 24 "Nesting" instabilities can 
occur for microscopically repulsive interparticle interac- 
tion strengths, and arise through the exchange of the 
nesting momenta {Kn} through the particle-hole chan- 
nel of these interactions. Such models may also exhibit 
the BCS superconducting instability,- 2 ^ which exists for 
any TRI Fermi liquid with (effectively) attractive pair- 
ing interactions. The ground state of such a half-filled, 
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FIG. 2: Brillouin zone (BZ) associated with the square lattice 
shown in Fig. [1] The shaded subregion is the sublattice BZ 
(sBZ) appropriate to the A and B sublattices. This subregion 
also indicates the Fermi sea at half filling, an attribute par- 
ticular to the square lattice model. The nesting wavevectors 
Ki and K2 are reciprocal lattice vectors for the sBZ. 
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Hubbard-like model with weak, but non-vanishing inter- 
actions is typically a Mott insulator, with lattice trans- 
lational symmetry spontaneously broken at the nesting 
wavelength, a superconductor, or a mixture of these, such 
as a "supersolid."- 2 ^ 2 ^^ 

A simple calculation 23 using the random phase approx- 
imation (RPA) predicts that the Fermi liquid phase of the 
(not disordered) model defined by Eq. |T]) is unstable to 
charge density wave (CDW) order for any V > U > 0. 
The CDW state is a Mott insulator, in which a greater 
proportion of the fermion density resides on one sublat- 
tice than the other. For the 2D case of the square lattice, 
the RPA calculation predicts a transition to the CDW 
state at a temperature 2 - 



2t exp 



t 



V-U 



(7) 



in the weak coupling limit < U < V <C t. Alternatively, 
a one-loop renormalization group (RG) calculation, 24 
performed upon a low-energy effective field theory de- 
scription of the model given by Eq. ([T|) , shows that the 
effective CDW interaction strength, which we define here 
as 



W C ~-(U-V), 
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(corresponding to a staggered charge density) grows to 
large negative values under renormalization if its initial 
value Wq " 1 < 0, i.e. was negative to begin with. This 
run off to strong interaction coupling is taken to signal 
the onset of CDW formation.™ These and analogous re- 
sults regarding the Neel ground state of the half-filled, 
spin-1/2 Hubbard model in d > 2 are well established, 
and the latter have been further confirmed with exten- 
sive numerical work. (See e.g. Refs. l22||27lf30l .) The spin- 
less Hubbard-like Hamiltonian in Eq. (TTJ), on the other 
hand, has received less attention in the numerical liter- 
ature; early Monte Carlo studies^ of the 2D version of 
this Hamiltonian, with U = 0, show the existence of a 
CDW ground state for positive V, but were unable to 
access the weak coupling limit V — ► + . 

Sublattice symmetry plays a crucial role in establish- 
ing the CDW ground state of the model in Eq. (fT]) (for 
V > U > 0). Nesting can also occur in the absence of 
SLS, since most aspects of Fermi surface geometry typ- 
ically depend strongly upon microscopic details. Under 
the repeated application of an appropriate renormaliza- 
tion group transformation, 24 however, this geometry is 
expected to deform continuously. SLS protects the nest- 
ing condition even as other details (such as the presence 
Fermi surface van Hove singularities) mutate through the 
RG process. The persistence of nesting and unbroken 
SLS, up to the onset of Mott or supersolid order, al- 
lows arbitrarily weak repulsive interparticle interactions 
to destabilize the Fermi liquid phase. 



B. Quenched disorder 

Now we turn to the incorporation of quenched (static) 
disorder into the model given by Eq. (TTJ. First, note 
that the addition of on-site (or: "diagonal") randomness, 
characterized by a strength Xp, breaks sublattice symme- 
try [Eq. ^] in every realization of disorder. Thus, turn- 
ing on diagonal disorder is expected to destroy the CDW 
ground state ; 31 ! 32 at least for sufficiently weak interac- 
tions. [Keeping the disorder strength Ad fixed, for ex- 
ample, we expect the absence of zero temperature CDW 
order in a window of (small) interaction strengths, e.g. 
for < V < V C (X D ), with {7 = 0)]. The proximity of 
a Mott insulating phase to the non-interacting (ballis- 
tic) Fermi liquid is the most essential characteristic of 
the Hubbard-like lattice model defined in Eq. ([1]). We 
conclude that this characteristic is lost upon the incor- 
poration of diagonal disorder. Studies of Hubbard-like 
models subject to diagonal disorder include those listed 
in Refs. l3l|]32ll33l and Refs. I3^l35l]36l for spinless and spin- 
1/2 lattice fermions, respectively. 

Instead, we consider the Hamiltonian in Eq. (fT]), 
weakly perturbed by purely "off-diagonal" randomness, 
which is taken to occur only in the intersublattice hop- 
ping amplitudes. Such a model is invariant under the 
sublattice symmetry (SLS) transformation [Eq. for 
each and every static realization of the disorder. We con- 
sider complex random (nearest-neighbor) hopping that 
breaks TRI. [Without interactions, the model is in the 
"crural" symmetry class AIII of Ref. 0-] The full lattice 
Hamiltonian is given by H = Ho + 5H, where Ho was 
defined in Eq. (fl]), and 
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Here, the random part of the hopping matrix element, 
St i j , is taken to be a Gaussian complex random variable 
with zero mean, independent on different lattice links. 
For our system of spinless fermions, this is consistent with 
the application of a random magnetic field to the other- 
wise clean model. For non-interacting spin-1/2 fermion 
systems, with nearest neighbor hopping on bipartite lat- 
tices at half filling, in random orbital and/or Zeeman 
magnetic fields, a general classification is (briefly) men- 
tioned in Sec. IIIB1 with more details being provided in 
Appendix [C] We will conclude from this classification 
that the results of this paper also apply to a related 
Hubbard-like model of spin-1/2 fermions, subject to an 
orbital magnetic field and to (homogeneous or random) 
spin-orbit coupling. Independent work on the effects of 
interparticle interactions upon these and several other 
classes of disordered spin-1/2 fermions has very recently 
appeared in Ref. H3- 

Our principal motivation for studying the model in 
Eqs. ((T]) and © is that, due to the presence of SLS, 
we expect both disorder and interparticle interactions to 
play important roles in the description of the low-energy 
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FIG. 3: Features of random hopping model (RHM) physics 
in 2D, in the absence of interparticle interactions. Subfigures 
(a) and (b) depict the qualitative energy (e)-dependence of 
the density of states v and inverse localization length l/£, re- 
spectively, in the chiral 7 - orthogonal (BDI) and unitary (AIII) 
class RHMs in d = 2. Both v and £ diverge upon approaching 
the band center, taken to occur at e = 0.— ' 12 ' 13 



FIG. 4: Same as Fig. [3] but for spatial dimensionalities d > 2, 
e.g. non-interacting random hopping model physics in 3D. 
Here, the shaded area represents a band of extended (delocal- 
ized) states; e c is the mobility edge. Increasing the strength of 
the random hopping disorder narrows the region of extended 
states in d > 2, as indicated by the arrows in subfigures (a) 
and (b). The DOS is finite, albeit parametrically enhanced 
at the band center.— 1 ^ 



physics. Because random hopping preserves the special 
SLS, our disordered model retains the "nesting" CDW 
instability of the associated clean system. This insta- 
bility can therefore compete with the unusual localiza- 
tion physics arising from SLS in the disordered, but non- 
interacting model (see below). The further assumption of 
broken TRI guarantees that we do not have to confront 
an additional superconducting instability) 3 - ' 22 ' 23 ' 24 ' 38 We 
note that the effects of hopping disorder upon the Neel 
ground state of the (slightly more complex) spin-1/2 
Hubbard model at half filling were studied numerically 
in Refs. and HO, although these studies were limited 
to d = 2. 

A second motivating factor is that, interestingly, and 
as alluded to above, the presence of SLS radically changes 
the localization physics of the disordered, non-interacting 
random hopping model [Eqs. (P) and ^ with V = U = 
0]. SLS enables the random hopping model (RHM) to 
evade the phenomenon of Anderson localization. Specif- 
ically, the non-interacting system exhibits a critical, de- 
localized phase at the band center (half filling) in one, 
two, and three dimensions for finite disorder strength, 
with a strongly divergent low-energy density of states in 
d = 1, 2 j s i 9 i 10 i 11 i 12 i 13 i 14 i 15 i 16 >11 In particular, there is no 
MIT and no Anderson insulating phase in d = 3 (in the 
absence of interactions). The essential features of random 
hopping model physics are summarized in Figs. [3] and 
|4j These figures apply, for example, to spinless lattice 
fermions with sublattice symmetry, with or without TRI 
(classes BDI and AIII, respectively) r- in d > 2. Upon ap- 
proaching the band center (e = 0) in 2D, both the density 



of states v(e) and the localization length ((e) strongly di- 
verge, as indicated in Fig. [3] By contrast, RHMs in d > 2 
support a diffusive metallic phase^S^ characterized by a 
band of delocalized states with energies |e| < e c , repre- 
sented by the shaded regions in Fig. HJ (e = ±£ c are 
particle and hole mobility edges.) While e c is expected 
to decrease monotonically with increasing disorder, as in- 
dicated by the arrows along the dashed boundaries of the 
shaded region in Figs.BJa) and[4jb), it is believed that in 
d > 2, there remains always a region of delocalized states 
of finite thickness in energy, centered around e = 0, for 
any finite disorder strength in a non-interacting system 
with SLS. 

Indeed, random hopping models have been of signif- 
icant theoretical interest in the recent past, both be- 
cause of the unusual delocalization physics described 
above, but also because these models have proven 
amenable to a variety of powerful analytical techniques 
in d < 2, with many exact and/or non-perturbative fea- 
tures now understood^ ' 13 ' 15 ' 17 This situation should be 
contrasted with our understanding of the conventional 
non-interacting ( "Wigner-Dyson" ) MIT, which is based 
largely upon perturbative analytical results in d > 2, us- 
ing the e-expansion4i 



C. Summary of results and outline 

In this work, we analyze the stability of the diffusive 
Fermi liquid phase of the Hubbard-like model defined by 
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Eqs. {1} and ([9]), in the simultaneous presence of both dis- 
order and interactions. We derive a low-energy, contin- 
uum field theory description of this system, which takes 
the form of a (class AIII) 7 Finkel'stein non- linear sigma 
model (FNLcrM)^ We employ the Schwinger-Keldysh 42 
method in order to ensemble average over realizations of 
the hopping disorder. The FNLerM contains parameters 
which specify the random hopping disorder and the inter- 
particle interaction strengths. We compute the one- loop 
renormalization group flow equations for these parame- 
ters, using a Wilsonian frequency-momentum shell back- 
ground field methodology.- We then discuss the physics 
of our model in (i) d = 2 and (ii) d > 2 dimensions. We 
now briefly summarize our results: 

For the 2D case, we find that the conducting phase of 
the disordered, non-interacting system is destabilized by 
the interactions. By contrast, such a phase, a metallic 
'diffusive Fermi liquid,' does exist (trivially) in d > 2; 
we identify what we call an "Anderson-Mott" disorder- 
driven instability of this metallic 'diffusive Fermi liquid.' 
This instability arises solely from the competition be- 
tween disorder and short-ranged interactions, and is per- 
turbatively controlled via an e-expansion in d = 2 + e. 
The instability that we find is indicative of a first-order 
metal- insulator transition (MIT); it occurs for disorder 
strengths of order e, and is therefore perturbatively acces- 
sible for e <C 1. We expect our result to be relevant for the 
Hubbard- like model in Eqs. (fTJ) and (O in three spatial 
dimensions (e = 1) for sufficiently strong disorder, and 
we stress that this "Anderson-Mott" instability is clearly 
distinct from the pure Mott "nesting" instability, which 
is driven solely by interactions (and not by disorder), al- 
though the latter also appears in the phase diagram of 
our model (see Fig. [231 below). The non-interacting ran- 
dom hopping model, Eqs. (|TJ) and © with V = U = 0, 
has no localized phase with disorder in the absence of in- 
teractions (see e.g. Section TV A ip ; therefore, a localized 
phase can only appear due to the presence of the interac- 
tions. The discovery of this disorder-driven, interaction- 
mediated diffusive Fermi liquid instability was previously 
announced in Ref. |5|. In this paper, we present a deriva- 
tion of this result, a detailed analysis of the phase dia- 
gram of the model, and we discuss our findings in view 
of previously known results for other, related disordered 
and interacting fermion systems. 

We note that the e-expansion is employed in this work 
as a technical tool in the continuum FNLerM descrip- 
tion, and should be thought of as a controlled approx- 
imation scheme to access the physics of the disordered, 
interacting Hubbard-likc Hamiltonian in Eqs. fl} and © 
in 3D. Clearly, we cannot easily define a bipartite lattice 
fermion model in a fractional number of d = 2 + e- dimen- 
sions; instead, we work with a continuum field theory, 
the FNLctM, argued to capture the low energy physics of 
the lattice model in both 2D (e = 0) and 3D (e = 1), 
whose internal structure is constrained by the crucial 
SLS [Eq. The field theory action of the continuum 
FNLerM [displayed in Eqs. ([57]) -(EH), Section ILTX4] can 



be analytically continued between integer dimensions in 
the usual way. The internal structure of the FNLctM, 
and thus SLS, is preserved under this continuation; it is 
SLS, then, that gives meaning to such an interpolation 
between lattice models in disparate integer dimensions. 
The "Anderson-Mott" instability of the diffusive Fermi 
liquid identified in this work occurs for perturbatively 
accessible, weak disorder strengths only for e < 1, i.e. 
in 2 < d <C 3. We conjecture that this instability also 
exists in the 3D FNLctM, and should therefore be found 
(i.e. via numerics) in the 3D Hubbard-like lattice fermion 
model. This conjecture cannot be directly proven here, 
since the instability of the 3D FNLerM, if it exists, would 
occur in the strong coupling regime. 

The organization of this paper is as follows: 
In Sec. HH we formulate a Schwinger-Keldysh— path 
integral representation for the lattice model given by 
Eqs. (J!) and ©; we then derive the low-energy, contin- 
uum FNLerM description. The technical content of our 
work appears in Sections IIIII and IIVI We set up our one- 
loop, frequency-momentum shell renormalization group 
calculation in Sec. IIIII specifying the parameterization 
of the FNLctM and stating the necessary diagrammatic 
Feynman rules. The actual one-loop calculation is chron- 
icled in Sec. IIVI We analyze and discuss our results in 
Sec. El The reader less interested in calculational details 
may skip Sees. IIIII and ITVl entirely, and proceed from the 
end of Sec. |TT]directly to Sec. El 

A variety of elaborations, extensions, and technical de- 
tails are relegated to the Appendices. In Appendix [A] 
we show how the structure of the continuum FNLctM for 
the model in the absence of interparticle interactions^ fol- 
lows immediately from a symmetry analysis of the non- 
interacting Keldysh action for the random hopping model 
[Eqs. HJ) and ©, with U = V = 0]. Appendix [B] details 
an expansion used in Sec. [Ill while Appendix [Cl describes 
the random matrix symmetry classification (along the 
lines of Ref. |43I ) of disordered, bipartite lattice models 
for spin-1/2 electrons. [See also Sec. Ill Bj . Appendix 
iDl provides a surprising alternative interpretation of the 
class AIII Finkel'stein NLctM, studied in this paper, in 
terms of the spin-1/2 quasiparticles of a spin-triplet, p- 
wave superconductor. This quasiparticle system may be 
defined directly in the continuum, without reference to 
a "microscopic" lattice model or an additional sublattice 
symmetry. Finally, Appendix [E] collects the loop inte- 
grals required in the RG calculation presented in Sec. IIVI 



II. FNLerM FORMULATION 

In this section we derive the class AIII Finkel'stein 
non-linear sigma model (FNLctM) description of the 
Hubbard-like lattice model given by Eqs. HI) and I©, 
using the Schwinger-Keldysh 42 method to perform the 
disorder averaging. The results of the derivation are pro- 
vided below in Eqs. ([571) . f58])- and and interpreted 
in the discussion following these equations. 
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A. Class AIII Finkel'stein NLctM 

1. Schwinger-Keldysh path integral 

To begin, we envisage a zero temperature, d + 1- 
dimensional real time-ordered (T-ordered) path integral 
Zt for the model defined in Eqs. (p} and ((9]). As usual, 
we need to normalize this path integral to unity in or- 
der to perform the ensemble average over realizations 
of the hopping disorder. We employ the Schwinger- 
Keldysh (Keldysh) methodj 42 which exploits the iden- 
tity 1/Zt = Zrp, where Zf is an anti-time ordered (T- 
ordered) path integral for the same model. We write the 
Keldysh generating function 42 
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where the non-interacting action is given by 
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with t 4 a = t + 5t 4 „-. The interactions reside in 
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The generating functional defined by Eq. (|10p is an in- 
tegral over the Grassmann fields c? A , c^ A , and c°g, c°g, 
defined on the A and B sublattices of the bipartite lat- 
tice. Note the locality of the non-interacting sector of 
the theory [Si in Eq. (jTTJ)] in frequency to, versus the lo- 
cality of the interacting sector [S2 in Eq. (|T2"]) ] in time t. 
In Eqs. HI]) and (TT2]), indices {i,i'} and {j,f} label A 
and _B sublattice sites, respectively, while the "Keldysh" 
species index a S 1,2 denotes the T-ordered (a = 1) and 
T-ordered (a = 2) branches of the theory. The number 
£ a takes the values 



density fields in Eq. (fT2"j) are defined via n Ai , B j{t) 
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The FNLcrM that we are after is a matrix field the- 
ory; it therefore makes good sense to introduce compact- 
ifying matrix notation at this early stage. We think of 
the field CAi/Bj {c-Ai/Bj) as a column (row) vector with 
a single "superindex," which is a direct product of fre- 
quency {u} and Keldysh {a} indices: e.g. cai — * c Ai (uj). 
In the Keldysh formalism, it is often useful to further 
divide frequency into a product of (discrete) sign and 
(continuous) modulus spaces: u> — \lo\ sgn(w). Next, 
we introduce two commuting sets of Pauli matrices: the 
matrix E m acts in the sgn(w) space, while the ma- 
trix acts in the Keldysh species (T/T) space, with 
m,n S {1,2,3}. Employing the conventional basis for 
all Pauli matrices, we identify E3 — > sgnw^.^f S a,a and 

6 -> £" <W< <5 a ' a ' , with <W = 2ttS(uj - wT). [e a was 
defined in Eq. (JT3J) .] Finally, we define the single particle 
energy matrix uj = \£i\ E3, with \uj\ —+ \uj\ <5 W , W < 5°" a ■ 

We make a change of variables c~Ai/Bj —> CAi/Bj£3, 
and rewrite the non-interacting sector [Eq. (fTT|) ] of the 
Keldysh action as 



Si = S? + SSl 
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The clean bipartite hopping model appears in the term 
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The momentum integration in Eq. (|15p is taken over 
the sublattice Brillouin zone (sBZ); e(k) is the clean en- 
ergy band structure. 45 Equation (fTS)) gives the most com- 
pact representation that we will use for the clean, non- 
interacting action. Here we have introduced a third set of 
Pauli matrices, <r m , me {1,2,3}, acting in the sublattice 
flavor {A, B} space, and we have also promoted momen- 
tum to an operator, k ► k. The column vector c car- 
ries indices in the momentum k, mod-energy \u>\, sgn(aj) 
(E), sublattice flavor (a), and Keldysh (£) spaces, i.e. 



"A/Bk 



[oj) with all indices displayed. Here, a = 1,2 
denotes the Keldysh species. The disorder is relegated to 
the perturbation 



SSi 



E 



CAi St z ■ C B j 



°Bj St* j CAi 



(17) 



1, a = 1 (T-ordered), 
-1, a = 2 (T-ordered). 



(13) 

2. Disorder averaging and Hubbard Stratonovich decoupling 



The factors i!; a r) sgn(w) in Eq. (JTTJ) are frequency in- 
tegration pole prescriptions appropriate to T- and T- 
ordered correlation functions, with 77 — > + . Finally, the 



We now ensemble average over realizations of the com- 
plex random hopping amplitudes {5t i ■} appearing in 
Eq. (|17p . In order to simplify the derivation of the 
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FIG. 5: Homogeneous pairing (dimer covering) of nearest 
neighbor sites on the square lattice. 



FNLcrM, we employ the following artifice: we assume 
that of the z c nearest neighbors bonds surrounding a 
given site i belonging to the A sublattice of the bipartite 
lattice under study, only one such bond is disordered. 
(z c is the coordination number.) We assume further that 
the same type of bond (specified by its orientation) is 
made random at each and every A site, thus allowing a 
unique, orientationally homogeneous pairing (dimer cov- 
ering) of the A and B sublattice sites of the bipartite 
lattice. An example of such a pairing is provided for 
the square lattice in Fig. [5] This seemingly pathologi- 
cal constraint upon the disorder distribution allows for 
the quickest derivation of the low-energy effective field 
theory4£ 

With such a pairing established between each A sub- 
lattice site i and its associated B sublattice site we 
have 



St 



0, 



(18) 



We take the "dimerized" bond amplitudes {St t to 
be Gaussian random variables, independent on different 
nearest-neighbor lattice links and identically dis- 

tributed with the following purely real mean and vari- 
ance: 



St, 



st d 



and 



St* ...M. 



(19) 



(20) 



respectively. The overbars in Eqs. (| 19|) and (|20[) denote 
disorder averaging. Although we will be ultimately in- 
terested in the limit of zero mean bond dimerization, it 
will prove convenient in the interim to retain 5td\ m 7^ 0. 
The superscript "(m)" appearing on the right-hand side 
of Eq. (|2"D|) stands for "microscopic," indicating that 
the quantity A (m) is defined at the lattice scale. Using 
Eqs. Q17 p -(|20 p . we define the disorder-averaged action 



SSi via 

e iSSl _ I | m j2 exp 



ieA 



with 



A i 



exp 



it 
j(*) 



jEKi(i)-^im| 2 +^ 1 
i 

+ iSS dim , (21) 



A, 




(Tl + 1(72 



A; 



(Tl - UT 2 



where we have defined 



Ai -> A 



a, a 
* j - a, a 



(22) 

(23a) 
(23b) 



In Eq. (|2ip . Tr denotes a matrix trace over indices in 
the sublattice flavor (a), frequency (u>), and Keldysh (£) 
spaces. The field Ai, introduced in Eq. (|2"Tj) . is a ma- 
trix of fermion bilinears; Eq. (j2"2")l details the sublattice 
flavor space decomposition of Ai in terms of the purely 
off-diagonal Pauli matrices <7\ and &2 ■ Although they oc- 
cupy the off-diagonal blocks of Ai, the fields Ai and A.,.-. 
actually describe fermion degrees of freedom residing en- 
tirely on the A and B sublattices, respectively, as shown 
in Eq. (f23|) . Ai and At,^ each carry frequency {u>,u>'} 

and Keldysh species {a, a'} indices. Finally, a non-zero 
mean bond dimerization <5idim [Eq. (|19[) ] gives rise to the 
following homogeneous term in Eq. (|2~lj) : 



5S d 



-Std\m 22 ( 



CAiCBj(i) + CBj(i)CAi) 



= -c\ai </>ij(k) + a 2 <Mk) 



(24) 



The second line of Eq. (I24| expresses SSdim in momentum 
space; here we have adopted the same compact notation 
employed in Eq. (|16p , wherein a\ and a? act in sublattice 
flavor space, and k is the (crystal) momentum operator. 
The functions </>i?(k) and 4>i(k) in Eq. (|24[) are real and 
imaginary components of the "dimerization" function 



iii(k) + #/(k) = -St^me* 



(25) 



In this equation, n is a unit vector pointing in the di- 
rection determined by the chosen bond dimerization [see 
Eqs. CBJ) and (US), and Fig. ©]. 

We have adopted a suggestive notation in Eq. |23|) 



to denote the pure sublattice fields Ai and A 



Con- 



sider the following spatially uniform deformation of these 
fermion bilinears: 



Ai -> U A AiU 



it 

JW 



^B^j(i)^A> 



(26) 



where U A and U B are independent unitary transforma- 
tions in \u)\ ® £ ® £ (frequency® Keldysh) space: 



C/t f/ A = f/^ B 



li 



(27) 
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The transformation in Eq. (|26|) is clearly a symmetry of 
the disorder-averaged action SSi, as can be seen from 
Eqs. (|2Tj) - (f24 |l . In fact, there is a direct relationship be- 
tween this transformation and the symmetry structure of 
the non- interacting Kcldysh action [Eqs. (|T6|) and (fT7|) ]. 
in every fixed realization of the static disorder; the con- 
nection is articulated in Appendix [A] Equation (|26|) sug- 
gests that we may regard as the "Hermitian ad- 
joint" of its associated nearest neighbor Ai, justified on 
length scales much larger than the lattice spacing. The 
identification At,^ ~ A\ implies the Hermiticity of the 

composite matrix field in Eq. (|23| : A\ ~ A4. 

Next, we decouple all four fermion terms appearing 
in the disorder-averaged and interacting sectors of the 
theory with bosonic Hubbard-Stratonovich fields. In the 
disorder-averaged sector [Eq. (f2"Tj) ] , we write 



(next-nearest neighbor) lattice sites {i,j} and 
vanishes otherwise. Now we decouple 



js 2 



where 



T>p exp 



v 0=1,2 J 



P a (t) 



dt 



^p^ {t )X^p*{t) 



p> a (t)N a (t) 



P a B {t) 



(34) 



(35) 



Again, we have suppressed all position space indices in 
Eqs. (|34 | and (f35 |) . 



exp 



A" 1 



EM 



VQ exp 



E Tr 



2A< m > 



as + q 1 a 1 



(28) 



The matrix field Qi is taken to be Hermitian, and purely 
off-diagonal in sublattice flavor space, i.e. 



with 







<J\ + 1(72 



t a, a' 



(30) 



[Compare Eqs. (JUJ and l[23j) ]. 

Turning to the interacting sector, we re-write Eq. Q12p 

as 

S 2 = - T / dtNTa i t )XN a (t), (31) 



where superscript "T" denotes the matrix transpose op- 
eration, and 



X 



Ux Vy 
Vy T Ux 



(32) 



is a matrix in sublattice flavor space (a), with elements 
involving the position space coupling functions x and y 
(defined below) , while the fermion sublattice densities are 
encoded in the vector 



N a (t) 



n%(t) 



(33) 



We have suppressed all position space indices in 
Eqs. (|3"T ]H [3"3" 1) . In Eq. ([32]), the function y -> y itj 
(x — > Xiji) equals unity for pairs of nearest neighbor 



3. Saddle point and gradient expansion 

Gathering together the homogeneous hopping 
[Eq. (fT6|) ]. mean bond dimerization [Eq. ([24])] . disorder- 
averaged [Eq. (|2"8")) ]. and interacting [Eq. (f3~4"]) ] pieces of 
our theory, we perform the Gaussian integral over the 
fermion fields, and arrive finally at the following effective 
field theory: 

Z= VQVpexv[-S Q -S p -SuEx]> (36) 



where 



ieA 



(37) 



and 



s P = E 

a=l,2 



Sdet = -Tr 



dtp Ta (t)X- 1 p a (t), (38) 



(ln^ + ^iQ + p]) 



(39) 



The trace in Eq. (|39() is performed over indices in the 
position, sublattice flavor (a), mod-energy \ui\, sgn(cj) 
(£), and Keldysh species (£) spaces. The operator Gq 1 
in Eq. (|3"9")) represents the inverse of the (Keldysh) single 
particle Green's function for the clean, non-interacting 
hopping model described by Eq. (|16[) . incorporating, in 
addition, the mean bond dimerization from Eq. 



G^ 1 = E 3 |w| + iritais - &i |e(k) + fo(k)J - & 2 ^(k). 

(40) 

The functions 4>R aid 4>i are real and imaginary parts of 
the "dimerization" function defined by Eq. (|25[) . 

We look for a spatially homogeneous saddle point so- 
lution to the action given by Eqs. (f3"T)) and (f3"9")) , in terms 
of the matrix field Q, with p = <pR — <pj — in Eqs. 
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and (|40|) (i.e. ignoring the interparticle interactions, and 
considering the limit of zero mean bond dimerization). In 
the low-energy /long-time limit \ui\ — > 0, the structure of 
the saddle point solution is determined by the pole pre- 
scription piece (the term proportional to rf) in Eq. (|40|) . 
We make the standard ansatz 



Qsp 



1 

27 



(41) 



with r the elastic scattering lifetime due to the disorder. 
Then the saddle point condition reduces to the usual self- 
consistent Born approximation (SCBA) for the elastic 
decay rate 1/r: 



1 



A" 1 



sBZ 



d d k 

(2-K) d 



£ 2 (k) 



1 



{2rf 



(42) 



For the 2D case of the square lattice (with its concomi- 
tant van Hove singularities at half filling) , Eq. (j42|) gives 



1 

27 



A" 1 



2nt 



In 



A (m) 



(43) 



in the weak disorder limit A <m) <C t 2 . 

Next we consider fluctuations about the saddle point 
solution given by Eq. (|41|L Dominant within the diffusive 
metallic phase are the long- wave length, low-energy Gold- 
stone (diffusion) modes that preserve the saddle point 
norm l/2r. These modes are generated by applying a 
slowly spatially varying generalization of the symmetry 
transformation in Eqs. ([26]) and (f27f to Qsp- (See also 
Appendix [Al) In terms of the sublattice flavor space de- 
composition [Eq. ([29]) ], we define 



Q(r) 



IT 



Q(r) 



4r 



(r), (44) 



where 



Q(r) = ^(r)Q SP ^(r), (45a) 
Qt(r) = L> B (r)Q S p^(r), (45b) 



and 



Alternatively, we note that Eqs. 1|27). (|45j) . and J46J) im- 
ply the unitary constraint 



Qsp = S 3 ^ 3 . 



(46) 



Qt( r )Q( r ) = i 



(47) 



Our FNLctM will therefore be a field theory of the unitary 
matrix Q(r) -> Q"£,(r). 

Finally, we assemble the action for the FNLctM from 
the components given by Eqs. ([37]) - ([39)) . The unitary 
constraint [Eq. (|4"T|) ] renders Sq in Eq. ([57)) an irrelevant 
constant, so we concentrate upon 5det in Eq. (|39p. Keep- 
ing only the most relevant terms in a gradient expansion, 



one obtains 



Sdet ~i2r J d d vTv [{%\Cj\ + i V ± 3 i 3 ) [Qt(r) + Q( r )]j 
+ (2v f t) 2 J d d rTr ( VQ f (r) • VQ(r) ) 



5 dim I q 
DET i" 



(48) 



where the terms S^'^Y and Sj are defined below in 

and 



Eqs. (gni) and (|5i|) . respectively. Eqs. 
(|5^|) are derived in Appendix [B] 

Versions of the "energy" (symmetry breaking) and 
"stiffness" terms generic to localization sigma models, 
but specialized here to the Keldysh class AIII FNLctM, 
appear on the first and second lines of Eq. (f48|) , respec- 
tively. Here, Vf is the average Fermi velocity (at half fill- 
ing), while the elastic scattering lifetime r is determined 
by the SCBA [Eq. @2J]. We emphasize that the struc- 
ture of these two terms, written out explicitly in Eq. (|48[) . 
is independent^ of the peculiar "dimerized" bond distri- 
bution [Eq. (fT8|) ] assumed for the random hopping. 

The first term on the third line of Eq. (|48p is a per- 
turbation arising from the presence of a non-zero mean 
bond dimerization 5tdim 7^ [Eq. ijTSJ)]: 



5 dim 
DET 



5U 



cTrTr 



Qt( r )VQ(r) 



(49) 



where the unit vector n, originally introduced in Eq. (|25|) . 
specifies the orientation of the mean bond dimerization, 
and c is a constant. Eq. (|49|) is derived in Appendix [Bl 
Consider the local vector operator 



iV$(r) = Tr Q+(r)VQ(r) 



(50) 



where $(r) is the U(l) phase of the unitary matrix field 
Q(r). Eq. (|49|) implies that iV$(r) provides a coarse- 
grained measure of the local bond dimerization, i.e. of 
the orientation of the strongest nearest neighbor hopping 
bonds within a neighborhood of size the elastic scatter- 
ing length, in the sublattice symmetric random hopping 
model. This interpretation can be argued on very gen- 
eral symmetry grounds, at least for the case of the hy- 
percubic lattice: here, a bond dimerization homogeneous 
in both magnitude and orientation preserves both SLS 
and sublattice translational symmetry, while breaking 
lattice rotation, reflection, and composite (intersublat- 
tice) translational symmetries: see e.g. Fig. [5j Eq. ([49]) 
is the most relevant term (in the sense of the renormal- 
ization group proximate to the diffusive metallic phase) 
consistent with these conditions. [A homogeneous back- 
ground dimerization, as in Eqs. ([24]) and (|25|) . breaks 
lattice reflection invariance through a lattice plane per- 
pendicular to the direction of the dimerization vector n; 
a composite sublattice translation involves an exchange 
of sublattice degrees of freedom, represented by Q <-> 
in the continuum theory.] 

Given this interpretation, we can generalize Eq. (|49|) 
to the case of a perturbation involving a static random 
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vector field n — » ndi m (r): 



where 



5'i 



dim 
DET 



cTrTr 



Q t (r)VQ(r) ■ n dim (r) 



(51) 



In Eq. (|5ip . ndi m (r) is taken to vary in both orien- 
tation and magnitude, and represents long- wavelength, 
quenched orientation fluctuations in the bond strength 
dimerization of the random hopping. Since Eq. (|5ip is 
consistent with the sublattice symmetry of the underly- 
ing lattice model, its effects should be included in the low 
energy effective theory upon the grounds of universality. 
We see that in order to obtain a final effective field theory 
that possesses rotational, translational, and reflection in- 
variances on average, it is necessary to specify a second 
measure of the random hopping strength, in addition to 
the parameter A (m) , the latter of which appears in the 
SCBA, Eq. dHJ), and was introduced in Eqs. O and 
dHJ). We take ridi m (r) to be a Gaussian random vari- 
able of vanishing mean, delta function correlated with 
the variance 



<Jr)4. m (r') = \ A 6 i >isW(r-r'). (52) 
Using Eqs. ([51]) and ([52")l . we define the disorder-averaged 



action S$fP T via^ 



_ edim 

e DET 



2?n dim exp 



2X A 



exp [*±J d d v [*V$(r)] 2 ] , 



(53) 



where the operator zV<I> was defined in Eq. f|50[) . In 
the low energy effective field theory, then, the random 
hopping is characterized by the two parameters A (m) 
[Eq. ((201)] and \ A [Eq. ([52"])]: the former sets the elastic 
scattering lifetime and the conductance ; 9 ' 12 ! 13 and there- 
fore reflects the "microscopic" structure of the random 
hopping up to small distance scales of order the Fermi 
wavelength (responsible for clastic backscattering events 
involving large crystal momentum transfers), while the 
latter characterizes the orientational fluctuations of the 
random hopping at larger distance scales (of order the 
mean free path). 

The last term on the third line of Eq. is due to 
the interactions, and is given by (see Appendix [Bj) 

S z = i 2r ]T fdt d d r [p\ (t, r) Q a $ (r) 

a=l,2 J 

+ ^(t,r)Qt r(r) ]. (54 ) 



Using Eq. (|38[) , we can now perform the Gaussian integral 
over the auxiliary field p, leading to the result 



Si -» i J2 ? I ,,ld " v 



a=l,2 



IV IV 



a, a • 
t,t 
f a, a 



(55) 



Tu ~ (2t) 2 U, T v ~ (2r) 2 U 



(56) 



[c.f. Eqs. pTjl and (f3"2"]) ]. and 1/r is the decay rate, 
Eqs. (g2J) and fig]). In moving between Eqs. ([3"5]). (|5l| . 
and (|55|) . we have approximated the short-ranged posi- 
tion space functions x and y, introduced as elements of 
the sublattice matrix X in Eq. (|32|) and defined below 
Eq. (J221), as x - y -> 8^ d \r - r'). 



4- The FNLaM and its coupling constants 

Combining the results of the previous subsection, 
Eqs. ([48]) . (|53f . and ([55| . we arrive at last to the final 
form of the Finkel'stein NLcrM (FNLcrM) description of 
the Hubbard- like model defined in Eqs. |T]) and (J9]). The 
FNLcrM is given by the functional integral 



Z = VQe 



-Sd — Si 



(57) 



where 



S ° = 1 kJ ( V ^ (r) ' V ^ (r) ) 

+ih J d d rTr ((S 3 p|+i7 7 S 3 e 3 )[Q t (r) + Q(r)]) 

/ d d r (Tr[Qt( r)V Q(r)]) 2 , (58) 



and 



* = l E ? I dtddr [2T V Q^' a (r)Q^ a (r) 

a=l,2 

+Y u [Q^{v)Q^{v)+Q\^{v)Q\^{v)]) . (59) 

The field variable Q(r) in Eqs. (|57| ([59")) is a complex, 
infinite-dimensional unitary matrix, 



Q T (r)Q(r) = 1, 



(60) 



carrying Keldysh species ({a, a'}) and time ({t, t'}) or 
frequency ({a;, a/}) indices: 

0(r)-Q:fJ,(r) or Q(r) - Q°;°'(r). (61) 

In Eq. ([57]) . 2?<3 is the invariant (Haar) functional mea- 
sure for the group manifold associated with Q. The ma- 
trix Q and its adjoint Q may be interpreted as contin- 
uum versions of the same- sublattice fermion bilinears 



(62a) 
(62b) 



which follows from Eqs. HHJ), (23J), (gSJ), and (25). It 
is a peculiar feature of the random hopping model that 



Q$ - c A {t)c A {t'), 

)l a / - c a B (t)c a B (f), 
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the low energy fluctuations on a given sublattice are de- 
scribed by a unitary, rather than a Hcrmitian matrix, 
familiar from other models. 

The Tr in Eq. ([58]) denotes a matrix trace over Keldysh 
species and time or frequency indices; the diagonal Pauli 
matrices S3 — > sgn uj 5^. u ' 5 a ' a and £3 — > £ a S^^> 8 a > a 
act in the sgn(w) and Keldysh species spaces, respec- 
tively, while \w\ — ► \ui\ S UtU > S a ' a is a matrix of abso- 
lute frequencies. Here and in Eq. ([55)) , the number £ a 
was defined in Eq. (|13p . The saddle point configuration 
Q(y) = QsP) determined above in Eq. ([46)) , is repeated 
here for convenience: 

Qsp = Ms- (63) 

The structure of Qsp is set by the symmetry breaking 
(Keldysh pole prescription) term, proportional to 77 — > 
0+, in Eq. ((55)) . 

The action given by Eq. ([55]) describes the low-energy 
diffusive physics of the non-interacting random hop- 
ping model^ ' 12 i 13 a replica version of the non-interacting 
sigma model with action Sd was originally studied by 
Gade and Wegner3 This non- interacting sector of the 
FNLcrM is parameterized by three coupling strengths: A, 
Xa, and h. The coupling constant A is a coarse-grained 
measure of the microscopic hopping disorder strength, 
A (m) [Eq. ([2TJ]) ]; we do not provide a precise relationship 
between A and A (m) , since such a connection depends 
upon various non-universal lattice-scale details. It can 
be shown that 1/A is proportional to the dimensionlcss 
dc conductance g of the systemi^i 3 -^ The parameter Xa 
gives a second measure of the disorder strength, unique 
to this "sublattice symmetry class," which strongly influ- 
ences the behavior of the low-energy, disorder-averaged 
single-particle density of states ^ 12 ' 13 [In Eq. {55]), we 
have rescaled A^ by a factor of 1/A 2 relative to Eq. ([53"]) . 
so that A and A ,4 now share the same naive "engineering" 
dimension. See also Sec. IIVBI ] The parameter A^ may 
be simply interpreted as characterizing the strength of 
long- wavelength, quenched orientational fluctuations of 
bond strength dimcrization in the microscopic random 
hopping disorder. [See the discussion following Eq. (14"5|) 
in the previous subsection for details.] Finally, the pa- 
rameter h in Eq. {55J is a dynamic scale factor, intro- 
duced here in order to track the scaling relationship be- 
tween length and time as the model is renormalizedi 2 ^ 3 - 

The interparticle interactions appear in the second 
term Si, given by Eq. ([55]). With the advent of Eq. ([52"]) , 
we may interpret <2"' t a (r) and Q\l' a (r) as continuum lo- 
cal density operators on the 'A' and l B' sublattices, re- 
spectively. Then the coarse-grained interaction strengths 
IV oc V and IV oc U in Eq. ([55)) couple to generic, 
short-ranged intersublattice and same-sublattice density- 
density interactions, respectively [compare to Eq. ([T]) 
above] . 

We use the renormalization group (RG) to study the 
model defined by Eqs. ([57)1 - ([55)) . It will prove convenient 



to define the linear combinations of IV and IV, 

r s = ±(iv + ±v), r c = ±(iv-rv). (64) 

In the continuum FNLcrM [Eq. ([55)) ]. the interaction 
strength T s couples to the squared (smooth) local charge 
density, [Q"; t a (r) + Q^' a (r)] 2 , while r c couples to the 
squared sublattice staggered charge density, [Q"' t °(r) — 

Q\ 1 ,a (r)] 2 . In accordance with the discussion in the para- 
graph below Eq. {5J in Sec. [I] we expect r c < to pro- 
mote charge density wave formation, while T c > should 
suppress it. 

We compute the one loop flow equations for the cou- 
pling constants A, A^, h, T s , and r c in the following two 
sections. We analyze and discuss our results in Sec. [V] 
The reader less interested in calculational details may 
skip Sees. IIIII and IIVI entirely, and proceed immediately 
to Sec. El 



B. Related models for spin- 1/2 fermions 

To conclude this section, we briefly discuss some con- 
nections between various models of disordered (and possi- 
bly interacting) spin-1/2 fermions, and their correspond- 
ing random matrix theory classification and sigma model 
descriptions. Details are provided in Appendices [U] and 

El 

Consider a clean system of spinless or spinful fermions 
with homogeneous, real nearest neighbor hopping on a 
bipartite lattice at half filling. In both the spinless and 
spinful cases, such a tight-binding model possesses three 
additional (non-spacetime) discrete symmetries: time- 
reversal invariance (TRI), sublattice symmetry (SLS), 
and particle-hole symmetry (PH). For spin-1/2 electrons, 
we also have spin SU(2) rotational symmetry. Here, we 
define the unitary particle-hole transformation as a prod- 
uct of antiunitary time-reversal and sublattice symmetry 
transformations; as a result, in the spin-1/2 case, the PH 
transformation involves a spin flip-see Appendix [C] for 
details. The introduction of quenched disorder may break 
or preserve each of these internal invariances, and the re- 
sulting disordered (non-interacting) Hamiltonian can be 
classified using random matrix theory.—^ 3 - 

Consider now the effect of random magnetic fields upon 
the otherwise clean spin-1/2 hopping model. Here, it is 
crucial to distinguish between the cases of random orbital 
and random Zeeman magnetic fields. A random orbital 
field preserves spin SU(2) rotational symmetry and SLS, 
while breaking TRI and particle-hole symmetry (PH) in 
every static disorder realization. A random Zeeman field, 
on the other hand, preserves PH, but breaks TRI, SLS, 
and spin SU(2) rotational symmetry (completely). These 
two cases actually fall into different symmetry classes 
of non-interacting, disordered quantum systems. In the 
classification scheme of Ref. 0, the random orbital field 
model, with SLS [and spin SU(2)] only, belongs to class 
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AIII, while the random Zeeman field model, with PH 
only, belongs to class C. These results are derived in Ap- 
pendix [C] A model with both random Zeeman and or- 
bital fields falls into the standard unitary Wigner-Dyson 
class A, since all three of the discrete symmetries TRI, 
SLS, and PH are broken. As discussed in the Introduc- 
tion, known (non-interacting) realizations of the "crural" 
class AIII^ i 10 ! 11 ! 12 ' 13 ! 14 ' 16 including the non-interacting 
version of the (spinless) model [Eqs. ([1]) and (|9|) with 
U = V = 0] studied in this paper, possess a conduct- 
ing phase with extended states in d = 2 spatial dimen- 
sions. A system of non-interacting, spin-1/2 supercon- 
ductor quasiparticles, with broken TRI, and with spin 
SU(2) rotational symmetry preserved in every instance of 
the static disorder, furnishes a better known realization of 
class Cj£ i 38 i 43 ' 49 i 50 i 51 i 52 Non-interacting systems belong- 
ing to classes C or A do not possess extended states in 
2D ; 41 i 49 i 53 except at their respective quantum Hal l 52 i 54 
transitions. 



From this discussion, we conclude that the FNLerM 
that we study in this paper also applies to a related 
Hubbard-like model for spin-1/2 fermions, subject to a 
random orbital magnetic field, and possessing strong, ho- 
mogeneous spin-orbit coupling. Spin-orbit coupling is 
needed to suppress an additional hydrodynamic spin dif- 
fusion channel, which we do not treat in this work. The 
2D half-filled, spin-1/2 Hubbard model subject to a ran- 
dom Zeeman field was studied numerically in Ref. |40| . 
The effects of interparticle interactions in the context of 
the superconductor quasiparticle interpretation of class 
C were analyzed using the FNLerM in Ref. [H; see 
also Ref. l5l|. Very recently, DellAnna^I has indepen- 
dently studied several universality classes of Finkel'stein 
NLerMs, including realizations of both the particle-hole 
symmetric class C and the sublattice symmetric class 
AIII for electrically charged spin-1/2 fermions; his cal- 
culations should capture the low energy physics of the 
above-described spinful Hubbard model in random Zee- 
man and orbital magnetic fields, respectively; in the lat- 
ter case, Dell'Anna includes the spin diffusion channel 
(i.e. assumes no spin-orbit coupling). His results are dis- 
cussed briefly in the conclusion to this paper, Sec. IV Dl 
Finally, we note that different classes appear in the pres- 
ence of TRI.ii 3 - 



Surprisingly, our AIII Finkel'stein NLctM equivalently 
describes a system of spin-1/2 superconductor quasipar- 
ticles, subject to disorder and interactions, with TRI and 
a U(l) remnant of the spin SU(2) rotational symmetry 
preserved in every static disorder realization. Such a 
(gapless) quasiparticle system could occur, e.g., in the 
polar phase of a p-wave, spin-triplet superconducto n 55 i 56 
This quasiparticle system may be defined directly in the 
continuum, without reference to a lattice model or an 
additional sublattice/chiral symmetry. The connection 
is derived in Appendix [D] 



III. PARAMETERIZATION AND FEYNMAN 
RULES 

We turn now to the setup of our perturbative, one- 
loop renormalization group calculation for the FNLerM 
defined by Eqs. (|57|) - (|59|) . The actual RG computation 
follows in Sec. lIVI Beyond developing the apparatus nec- 
essary for the RG, the material in this section serves also 
to further elucidate the structure of the sigma model de- 
scription of sublattice symmetric disorder and interpar- 
ticle interactions. 

To begin, we shift the saddle point £3^3 [Eq. (|6"3")) ] to 
the identity 1 via left group translation of the unitary 
matrix field Q(r), 

Q(r)^± 3 £ 3 Q(r), Qt( r ) _> Q\r)£ 3 £ 3 . (65) 

The FNLerM action, Eqs. iJSSJ and (|55|). becomes 



vg f • vg 



Sd =2X ' d rTl 



^/A [Tr(gtvg 



ihf d d rTi [(\Cj\i 3 + iiv)(Q 1 +Q)) , (66) 



and 



dui du>2 du>3 du>4 w] . 

o r 01+3 9+4 

2tt 2tt 2tt 2tt 1 + j ^+ 4 



(' 



■jf a,a^-\f a, a 
53.4 



tr s 1 s 3 Q 12 Q 34 + s 2 S4Q[ i2 g 3 
2r yS2S3 Q t li rQ^j. (67) 



Eq. (f6"T|) expresses Si in frequency space, where we have 
implemented the compact notations 



si = (%)ui,ui = sgn(wi), 



(68) 



etc., i.e. numeric subscripts represent associated fre- 
quency labels. 

The saddle point shift defined by Eq. (|6"5)) modifies 
only the symmetry-breaking term, proportional to h, in 
the non-interacting sector of the theory [c.f. Eqs. ([58]) 
and (|6"6")) ]. By contrast, the transformation in Eq. (|63|) 
inserts explicit factors of SiSj = sgn(wi) sgn(wj), i,j € 
{1, 2, 3, 4}, into all terms inhabiting the interacting sector 
Si [Eqs. (|59|) and (|67|) ]. We show below that these factors 
function as projection matrices, dividing the interaction- 
dressed diffusion modes into smooth and sublattice stag- 
gered charge density diffuson channels, characterized 
by interaction strengths T s and r c , respectively. [See 
Eq. dMD, as well as Eqs. fgQ ]) -(P ]) . below.] 

We now parameterize the FNLerM for the renor- 
malization group (RG) computation. We employ a 
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the origin \uj\ 



Dk 2 = 0. Here 



A A 



FIG. 6: Frequency-momentum shell of thickness (A — A) sup- 
porting fast mode field coordinates Y°[ ,a , (k). The cubic vol- 
ume of linear size A enclosed by the shell supports the slow 



mode fields <5Qc ' a 



,(k) and SQl'j^Ot). 



Wilsonian frequency-momentum shell, background field 
methodology^ The first step is to split Q into "fast" Q F 
and "slow" Q s mode parts, 



Q(r) = Q F (r)Q s (r) 



= Q f (t) 



5 (r) 



(69) 



where both fast and slow matrix fields satisfy the unitary 
constraint [Eq. (J5DJ)]: 



(70) 



On the second line of Eq. (|69|) . we have further 
decomposed^ the slow mode field Q s into the homo- 
geneous saddle point 1 plus the "small" perturbation 
SQ s (r) ee Q s (r) - i. 

We require a parameterization for the unitary fast 
mode field Qp in terms of some unconstrained coordi- 
nates; we choose geodetic coordinates on the group man- 
ifold, 



Q F (r) 



exp 



iY(r) fa 1 + iY{r) + . . 



(71) 



where Y> = Y is a Hermitian matrix belonging to the 
Lie Algebra that generates Q F . 

We explain now the meaning of the fast-slow decompo- 
sition defined by Eq. (f6"9")> . The spatial Fourier transform 
of the slow mode fluctuation 

SQ s (k) = [ d d r cx P Hk.r)<5Q s (r) - SQ^' ,(k) (72) 



is taken to possess support within a cube of linear size A, 
in the (3 dimensional) absolute frequency and squared- 
momentum space |o/|,-Dk 2 ); the cube encompasses 



D = l/(\h). (73) 
is the effective jiiffusion constant. The linear size of the 
cube satisfies A -C 1/r ~ D/l 2 (measured in energy 
units), where r is the elastic scattering lifetime and I the 
mean free path. We take the support of the fast mode 
coordinate 



Y (k) = J AexpHk • r) Y(r) - Y^,(k) 



(74) 



to lie within a thin frequency-momentum shell enclosing 
the slow mode cube, 

A < Dk 2 < A, < \lu\ < A, < < A; (75a) 
A < \lo\ < A, < Dk 2 < A, < < A; (75b) 
A < < A, < Dk 2 < A, < \lu\ < A, (75c) 

where A/A pa 1 + 2dl is a ratio of energy cutoffs, with 
< dl <C 1. The regions of slow and fast mode support 
are illustrated in Fig. [SJ In the next section, we will 
integrate out the degrees of freedom in this shell and 
determine the resulting effect upon the slow modes SQ S . 

Substituting the fast-slow decomposition [Eq. ([69]) ] 
into the action in Eqs. (|66|) and (|67|) . and expanding to 
second order— in the fast mode coordinates Y(r) using 
Eq. ([7T]), the Keldysh generating functional Z [Eq. ([57]) ] 
may be re-written as 



Z~y VQsVYcxp ^-S[Y,5Qs\) 



(76) 



We divide the action 5* in Eq. ([76")) into terms containing 
(i) only slow modes, (ii) only fast modes, and (iii) both 
fast and slow modes: 

S = S S [SQ S ] + S F [Y] + S F/S [Y,6Q S }. (77) 

The pure slow mode sector of the theory is 

S S [SQ S ] = S D [Q^1 + SQ S ] + S T [Q -» 1 + 8Q S ], (78) 

with Sd and Si as given by Eqs. (|6"6"j) and (fBTj) . Of the 
remaining terms in Eq. (|77[) . Spfy - ] contains only fast 
mode degrees of freedom, while Sf/s{Y, SQ S ] describes 
the coupling between the fast and slow mode fields. 

To one loop order, we need only retain terms in 5V [5^] 
[Eq. ([77)1 ] to second order in Y(r), thereby obtaining the 
fast mode Gaussian diffusion propagator. We write 



S F EE Sf > 



(I) 
F ' 



(79) 



where the sector of the fast mode theory independent of 
the interparticle interactions is 
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4 D) =-. 



Tr 



3 + liv)Y 



X 2 



Tr (vf) 



while the interactions give rise to the term 

^ = -E^ / ppppd^v (5 1+ 3, 2+ 4[r s ( Sl - S2 )( S 3- S 4) + r c ( Sl + S2 )( S 3 + S 4)]F-y 3 7) 

z — ' J 2,71 2,71 2,71 2,71 \ ' ' ) 



(80a) 



(80b) 



We have expressed S\ , Eq. (|80b[) containing the inter- 
particle interactions, in terms of the smooth and of the 
sublattice staggered (CDW) charge density interaction 
parameters T s and T c , defined by Eq. In Eq. (|80b|) . 

r s and r c couple to the frequency index tensors (s% — 
S2)(s 3 - s 4 )S 1+3t2 +4 and (si + s 2 )(s 3 + s 4 )<5i +3i 2 + 4, re- 
spectively, which project out the channels of the prop- 
agator (Y\2^t) off-diagonal and diagonal in sgn(w) 
space, respectively. These channels are orthogonal, so 
that Eqs. (j80a|) and (|80b|) may be simply inverted to ob- 
tain the fast mode propagator for the theory. Note that 



all Keldysh indices must be identical to obtain a non-zero 
contribution from Eq. (|80b[) . 

The fast mode propagator cleanly decomposes into four 
disparate components; in terms of the spatial Fourier 
transform Y(k) defined via Eq. ([74]) . 

(Y$(-k)Y 3 c : *(k)) ee P A + P A + P S + P C , (81) 

where the four symbols Pa, Pa, Ps, and Pc represent the 
components 



P A = A5j i (|wi|,H,k) S a ' d S c ' b 6 1A S 3 , 



Pa = - 



^k'Ag-dwiUwxI.kJA^dwal.lwal.k) 



6 a ' b 5 c > d 5i 2 (5; 



3.4 



Ps = 
Pc = 



2i£T s (si - s 2 )(s 3 ~ 84.) A%{\oj 1 - w 2 |, k)A^ Q (| Wl - w 2 |,0, k) 

2l£ 1 C (S1 + S2)(S 3 + S4) 



S^ S a,d S c,d Si+32+4 



l + lc f a {\u Jl -wa|,k) 



5 O '^ ' d 5 C ' rf 5l+3,2+4. 



(82a) 
(82b) 

(82c) 

(82d) 



Keldysh {a, 6, . . .} and frequency {1, 2, . . .} indices in 
Eq. ([52")) should be matched to those in Eq. (|8"Tj) . Com- 
ponents P A and P A [Eqs. (fg2"aj) and (|82b|) ] follow from the 
inversion of the disorder-only action in Eq. (|80a[) , while 
Ps and Pc [Eqs. (|82cj) and J82d|] incorporate the inter- 
actions from Eq. (|80b|) . 

The basic (heat) diffuson kernel in Eq. (|82ap is 

A£ fc (M, Iwal.k) = \ [D\?-i{e\ui\ + e\"2\)T 1 , 

(83) 

with the diffusion constant D defined via Eq. (175)) . The 
diffuson matrix propagator Pa in Eqs. (|82a[) is pictured 
as a thick line segment in Fig. [7(a), with ends that split 
into pairs of directed thin lines. As shown in Figs.UJd)- 
[7(f), each such thick line segment in Figs. E(a)-[7(c) can 
be taken to represent a pair of single fermion particle and 
hole lines, carrying counterpropagating arrows indicating 
the flow of the conserved electric U(l) current [See, e.g., 
Eq. (|62|) ]. The numeric labels appearing along the termi- 



nating thin lines in Fig. [TJa) encode the frequency and 
Keldysh indices carried by the fast mode matrix fields 
Y — ► Y*' u i . These indices are propagated along (unbro- 
ken) thick line segments without mixing, as in Fig. [7(d). 
In both Figs. [7] and [51 we employ the following conven- 
tion: different numerical labels encode independent fre- 
quency indices, and numerical labels carrying different 
numbers of primes indicate independent Keldysh species 
indices; numerical labels with the same number of primes 
represent indices that share the same Keldysh species. 
Each thick line splitting in Fig. [7J corresponds to the in- 
sertion of a fast mode field Y a ' a ,, , where the "left" indices 
{ui, a} accompany the arrow flowing out of the thick line, 
while the "right" indices {a/, a'} accompany the arrow 
flowing into the thick line. 

Inversion of Eq. (|80ap also gives the component Pa 
proportional to Xa, Eq. (|82b[) . which has a different fre- 
quency and Keldysh index structure than the basic dif- 
fuson P A [Eq. (pZaj) ]. P A is depicted in Fig. |7jb). The 
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function 



FIG. 7: Feynman rules I: fast mode propagator, from 
Eqs. (|81[1 and (|82[1 . Components of the disorder-only sector 
P\ and Pa are depicted in (a) and (b), respectively; (c) rep- 
resents the sum of interaction-dressed components Ps + Pc, 
which vanishes in the non-interacting limit. Subfigures (d), 
(e), and (f) show that each thick line in (a), (b), and (c), re- 
spectively, can be understood as a pair of particle and hole 
lines, which carry counterpropagating arrows to indicate the 
flow of the conserved electric U(l) current. Numeric labels 
represent the frequency indices of single fermion lines, while 
primes denote the associated Keldysh species indices. Nu- 
meric labels with the same number of primes share the same 
Keldysh index, while labels with different numbers of primes 
possess independent Keldysh indices. 



inversion of Eq. ()80a[) gives only the sum of Pa and Pa, 
i.e. terms up to first order in A^. We can try to build 
terms higher order in Xa by cascading together multi- 
ple such sections, but it is clear from Figs.[JJb) and[JJe) 
that such a construction necessarily contains at least one 
closed "Keldysh" loop, defined here as a simultaneous 
Keldysh species index summation and frequency integra- 
tion along a closed single fermion line. As with the replica 
trick, such closed loops vanish in the Keldysh formalism. 

The smooth and sublattice staggered interparticle in- 
teractions in Eq. (|80bp . characterized by T s and r c , 
respectively, dress the bare diffuson propagator Pa 
[Eq. (|82ap j. giving rise to the propagator components P$ 
and Pc, Eqs. (|82c[) and (|82d[) , respectively. These equa- 
tions show that the component Ps (Pc) projects upon 
the propagator channel off-diagonal (diagonal) in sgn(w) 
space. The kernel 

A a s (\u J \,k)^j i [Dk 2 - l (l- ls )a^\r 1 (84) 

appears in Ps, the channel of the interaction-dressed 
propagator off-diagonal in sgn(w) space [Eq. (|82cp ]. and 
is related 2 ^ to the diffusion of the physical, conserved 
electric U(l) charge. In Eqs. (f84|) and (182d|) . we have 
introduced the relative interaction constants 



4 


2 


irh 6 


irh 


_ 4 


2 


irh c 


Tth 



7c 



where we have used Eq. 



r, 



rv), 



(85a) 
(85b) 



/ a (M,k)^ln 



2A 



i£ a Dk z 



(86) 



Finally, the logarithmic 



with A a hard frequency cutoff [c.f. Eq. ([75)) ]. appears 
in Pc, the channel of the interaction-dressed propagator 
diagonal in sgn(w) space [Eq. (|82d[) ]. The structur o 3 ' 59 i 60 
of Pc does not follow from any conservation law, as the 
sublattice staggered electric charge density (the 'CDW 
order parameter associated with the interaction strength 
F c ) does not represent a conserved quantity. In this pa- 
per, we will work only to the lowest non-trivial order in 
7 C ; practically, this means ignoring 2 ^ the logarithmic 
denominator [Eq. ||SBJ|] in P c [Eq. (|82d])]. This approx- 
imation is adequate for all of the results presented in 
Sec. El The sum of the interaction-dressed propagator 
components Ps + Pc will be depicted as a thick line with 
a black dot, as shown in Fig. \T[c). 

The channel of the fast mode propagator diagonal 
in sgn(w) space, including the components Pa, pro- 
portional to the disorder strength Xa [Eq. (|82b[) ]. and 
Pc, proportional to the CDW interaction strength r c 
[Eq. (|82d[) ]. are special to a system with sublattice sym- 
metry [Eq. (HI)]. (Recall from Sees . Ill A 31 and III A 41 that 
Xa measures the strength of quenched orientational fluc- 
tuations of bond dimerization in the intersublattice hop- 
ping disorder.) The addition of SLS breaking, e.g. in 
the form of on-site disorder, changes the random matrix 
class^ of the (non-interacting version of the) random hop- 
ping model studied in this paper, Eqs. (JTJ) and © with 
U = V = 0, from the "sublattice/chiral" class AIII to the 
ordinary unitary metal class A. Equivalently, breaking 
SLS reduces the size of the sigma model target manifold, 
as discussed in Appendix [X] The crossover is marked by 
the appearance of a "mass" in the FNLtrM propagator 
channel diagonal in sgn(w), which gaps out the propaga- 
tor components Pa and Pc [as well as "half" of the basic 
diffuson modes Pa in Eq. (|82ap ]. This is the reason that 
the coupling strengths Xa and T c do not appear in the 
spinless unitary class Finkel'stein NLcrM^ which retains 
the parameters A, h, and T s . 

The fast and slow mode fields are coupled together by 
the term S F/S [Y, 5Q S ] in Eq. ([77]). S f/S [Y, SQ S ] gives 
rise to the Feynman vertices pictured in Figs. [51 We 
will refer to the vertex pictured in Fig. [8j60 as "V(6>)," 
with 9 E {a,b,c, ...} In these figures, the vertex cor- 
ners adorned by thick line "spokes" that split into di- 
rected thin lines represent insertions of the slow (SQ S 
and/or SQ^) or fast (Y) mode matrix fields. Specifically, 
the triangular black and blank "spoke" terminals denote 
the slow mode fields 5Q S and SQ^, respectively, while 
the half-black, half-blank terminals denote linear combi- 
nations of these. Each unterminated thick line "spoke" 
in Fig. [8J represents a fast mode field Y, which may be 
interconnected in pairs using the propagators shown in 
Fig. [7J with amplitudes given by the associated expres- 
sions in Eq. (152"[) . As in Fig. [7J we use numerical labels 
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FIG. 8: Feynman rules II: vertices coupling together fast and 
slow modes arising from the term Sf/s\Y , 5Q S ] in Eq. (|77|l . 
Vertex V(#) is pictured in subfigure (9), with 6 6 {a, b, c, . . .}. 
In this figure, black and blank triangular terminals indi- 
cate 8Q a and 8Q\ slow mode fields, respectively, while half- 
black, half-blank terminals indicate linear combinations of 
these. Unterminated thick lines represent fast mode fields 
Y . Vertices V(a) and V(b) derive from the stiffness and en- 
ergy terms, respectively, of the non-interacting sigma model 
action [Eq. l|66jl ]. while vertices V(c)-V(k) arise from the in- 
terparticle interactions [Eq. (|67[) ]. Momentum labels q, k, 
and q + k have been furnished for the stiffness vertex V(a), 
as this vertex depends explicitly upon these momenta. 



to indicate frequency indices, with primes to distinguish 
independent Keldysh species indices. At a given fast or 
slow mode field (thick line "spoke") insertion, e.g. the 
6Q$ — > <5Qs5 4 s l° w mode associated with the black tri- 
angular terminal in Fig.[5Jc), the "left" indices {5, a'} ac- 
company the arrow flowing out of the thick line "spoke," 



into the vertex, while the "right" indices {4, a} accom- 
pany the arrow flowing into the thick line, out of the 
vertex. 

The factors associated with the vertices 
{V(a), V(b), . . .} are listed in Table [fl In this ta- 
ble, the fast mode fields Y have been amputated from 
the vertex expressions; the structure of the fast-slow 
mode coupling should be understood from Fig. (|8j). 
Vertices V(a) and V(b) obtain from the non-interacting 
sector of the FNLcrM action, Eq. (|66| . while vertices 
V(c)-V(k) derive from the interparticle interactions, 
Eq. ([17]), indicated by the black dots in Figs. [U^-^k). 
The first entry of Table U gives the factor associated with 
the "stiffness vertex" V(a), shown in Fig. [5Ja). Here we 
have introduced the slow mode vector operator 



L s (r) = Q s (r)V^(r) = [1 + 5Q s (r)j VSQl(r). (87) 

Note that all fast-slow vertices pictured in Fig. [5J ex- 
cept V(f), are bilinear in the fast mode fields V(r) (rep- 



it, 



TABLE I: Feynman rules III: factors associated with vertices 
{V(a), V(b), . . .}, pictured in Fig. [5] coupling together fast 
and slow modes. Frequency {1, 2, 3, ... } and Keldysh species 
{a, a', a"} indices in the factor V(#) should be matched to 
those in the corresponding Figure [$!#)• Fast mode fields Y 
have been amputated from the vertex factors in this table 
for brevity; the fast-slow mode coupling structure may be 
understood from Fig. [8] The (spatial Fourier transform of 
the) slow mode operator Lj' 1 a 2 (q) is defined by Eq. (|87[1 . 

V(a)=-(i/2A)/ [d d q/(2^) d ] (2k + q) -L|-' 2 (q) 

V(b)=(ih/2) / d d r (\u»\e'6Qg a + \m\eSQiV,2) 

V(c) =2^1+3,2+4 (r uSl s 3 -FvS2S 3 ) J d d r5Q^ 4 

V(d) =2^1+3,2+4 {Y uS 2Si - V V S2S3) J d d T 5Ql ^ 

V(e) =iC5i+3,2+± I d d r [r,(*i + s 2 )(s 3 8Q^ A + s 4 5Ql^) 

+r c (si - S2)(s35Qs3% ~ S 4 5( 3s3,4)] 

=2e Q 5i+3,2+4/d d r [r s (si - s 2 )(s 3 6Q%f A + s^SQl^) 
+r c ( Sl + s 2 )(s 3 5Q°>« 4 - s 4 SQlt:t)] 
=2^ a 5i + 3,2 + 4rc s lS3 / d d r 5Q a s '<%5Q a s ".<2 
=2it a 6i +3 ,2+4,TuS2S4 J d d r SQl 1% &Q\ $$" 



V(f): 

v(g) : 

V(h): 
V(i) : 

v(i) = 

V(k) = 



-2iCSi+3a+iFvS2S 3 f d d r8Ql^' 5Qse]2 
^i+3,2+4 / d d r [rusis 3 5Q^%6Q^ 4 

r^i+3,2+4 / d d r [ruS2S45Ql a ^'SQl a 3 l2 

+V V S2S 3 SQll^ 5Qs3%] 
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resented by unterminated thick lines); the latter vertex 
involves only a single fast mode field. Such a vertex 
would vanish in a pure momentum shell treatment, in 
which one integrates undetermined fast mode loop mo- 
menta over a thin shell, while simultaneously integrating 
undetermined fast mode loop frequencies over the entire 
real line. We show in Sec. IIVI that a diagram involv- 
ing two copies of V(f) produces the term in the RG flow 
equations that gives rise to the CDW instability in the 
clean Hubbard-like model [Eq. {I}]. As the CDW term 
must be present in the advent of sublattice symmetry 
[e.g. in the ballistic limit A, Xa — > 0], we are forced to 
work with the frequency-momentum shell method in this 
paper. A similar term responsible for the BCS supercon- 
ducting instability of the diffusive Fermi liquid also arises 
in the FNLctM description of normal, TRI metals only in 
the frequency-momentum shell scheme.— ^ 59 ' 60 



IV. ONE-LOOP CALCULATION 

With the setup outlined in Sec. IIIII complete, we com- 
mence here the renormalization group (RG) calculation 
proper for the FNLtrM originally defined by Eqs. Q57[l- 
(|59[) . In two dimensions, the disorder parameters A and 
Xa, as well as the ratios 7 S and j c , defined by Eq. (|85|) . 
carry zero "engineering" dimension (demonstrated ex- 
plicitly in Sec. HVBl below). In the language of the RG, 
we say that these parameters are marginal (in d — 2) at 
tree level. In order to understand the infrared physics 
of the theory, we must therefore go beyond the reach of 
dimensional analysis. We compute here the one-loop RG 
flow equations for the coupling strengths A, Xa, and h, as 
well as Tu and Ty, or equivalently, T s and r c [Eq. 
The loop expansion is performed in d = 2 + e dimen- 
sions, with < e -C 1, and is formally organized as an 
expansion in powers of the disorder strength A, which 
is inversely proportional to the dimensionless DC con- 
ductance. As discussed below Eq. ([56"]) . we choose to 
work only to the lowest non-trivial order in the CDW 
interaction parameter r c (or 7 C ); by contrast, the loop 
expansion incorporates contributions from Xa, h, and T s 
to all orders. We need make no assumptions about the 
smallness of these latter three parameters. The diagram- 
matics appear in Sec. IIV At the results of this section are 
combined with dimensional analysis in Sec. lIV 51 yielding 
the desired one-loop flow equations. The flow equations 
obtained here are summarized and interpreted in Sec. [V] 



A. Renormalization 

The calculation is performed in the frequency- 
momentum shell background field formalism established 
in the previous section. We integrate out the fast modes 
Y££',(k) [Eq. (JTIJ] lying within the shell defined by 



Eq. ((751), where 

t re 1 + 2dl, (88) 
A 

and < dl <C 1 , as below Eq. ([75)) . [The quantity dl may 
be understood as an infinitesimal change in the log of the 
length scale L, with L the linear system size. In the dif- 
fusive metallic regime, energy lu scales like the square of 
an inverse length: u> ~ D/L 2 ; hence, the factor of "2" in 
the definition ([55)) .] The fast mode integration produces 
corrections to the action Ss [Eq. ([78)1 ] for the slow mode 
fields 6Qs^,(k) and SQ^,(k) [Eq ©]. The regions 
of fast and slow mode support were shown in Fig. [5] Note 
that in the case of the fast mode field Y"'"'(k), at least 
one of the three variables a € {Dk 2 , \cj\, is always 
"fast," lying within the range A < a < A. 

In the pure slow mode sector of the theory, with ac- 
tion Ss defined via Eqs. ([78]). (|66|) , and (|67|). all five pa- 
rameters A, A^, h, IV, and Ty couple to local opera- 
tors which contain non- vanishing terms quadratic in SQ$ 
and/or 5Q 5 £^ In order to renormalize the parameters of 
the theory, then, we need only consider corrections to 
terms up to second order in SQ S (and its adjoint) in the 
pure slow mode FNLctM action. 

The vertices necessary for the one loop calculation are 
provided by Fig. [5J and Table [U To compute all one loop 
corrections up to second homogeneous order in SQ$(r) 

and <5Qg(r), we must calculate diagrams involving one 
copy of each of the vertices pictured in Figs. [5Ja)-[5Jk), 
diagrams involving pairs of the vertices (in two graph 
topologies) in Figs. [8j(a)-[8je), and finally a diagram in- 
volving two copies of the vertex depicted in Fig. [5Jf). 
Taking into account the three propagator sections shown 
in Fig. [7J a naive estimate gives 3 x 11 + 3 x 3 x 15 x 2+3 = 
306 possible diagrams! However, many of these turn out 
to individually vanish, give corrections to higher order 
in r c , or correct only irrelevant operators in the FNLctM 
action. Among the diagrams that individually vanish are 
those that contain one or more closed "Keldysh" loops. 
A Keldysh loop is defined as a simultaneous Keldysh 
species index summation and frequency integration along 
a closed single fermion (thin) line. The contributions 
from the two Keldysh species exactly cancel in such a 
loop, as required by the normalization condition Z = 1. 
[See also Eq. (fT0]l and the discussion preceding it.] 

In this subsection we furnish the subset of 83 non- 
vanishing diagrams that correct the marginal parame- 
ters in the FNLctM and/or correspond to the genera- 
tion of marginal operators not originally present in the 
FNLctM action. All contributions of the latter type must 
cancel if the theory is to be renormalizable to one loop; 
we will show that this is indeed the case. The required 
83 diagrams are organized into 14 categories, shown in 
Figs. fTDI [22] renormalizing each of the coupling strengths 
appearing in the disorder-only and interacting sectors of 
the FNLcrM. To each category of diagrams, we provide 
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the label 23m, with m 6 {1, ...,14}. Category labels 
{231,232,...} appear in the captions of the associated 
figures {[TOlfTTl . . . }. Individual diagrams and their asso- 
ciated amplitudes will be referred to by a category label, 
a letter, and if necessary, a subscript Roman numeral. 
We give two examples: 231(a) refers to the diagram la- 
beled (a) in the left-hand column of Fig. [TO] (category 
5)1), while 2)2(b)j;j refers to the diagram labeled (iii) in 
the right-hand column of Fig. [IT] (category 232). Below 
we examine each category of corrections in turn. We 
also calculate the one loop renormalization of the sin- 
gle particle density of states (DOS) v(ui). Given the 
rather large number of non- vanishing diagrams, we ex- 
plain in detail the computation of only a handful of the 
associated one loop corrections. The goal here is to il- 
lustrate the process; calculation of the remaining dia- 
grams is straight-forward, if time consuming. Explicit 
frequency and Keldysh indices (using the same conven- 
tions employed in Figs. [JJ and [5]) distinguish the diagrams 
in Figs. [TOH2"2l whose detailed evaluation is provided in 
this section. 

The key ingredients for the renormalization process are 
the fast mode propagators Pa, Pa, Ps, and Pc, depicted 
in Figs.[JJ[a)-(c), with the associated amplitudes provided 
by Eq. (|8"2")) . and the vertices V(#) coupling together fast 
and slow modes, with 9 £ {a, b, c, . . .}, cataloged in Table 
U and pictured in Fig. [5J Explicit formulae for the nec- 
essary frequency-momentum shell loop integrations are 
relegated to Appendix [El 



1. Propagator with a twist 

Before we begin, we need to introduce one additional 
piece of diagrammatic notation. The basic diffuson prop- 
agator P\, given by Eq. (|82a[) . is represented as the thick 
line segment shown in Fig. [TJa). As shown in Fig. [TJd) 
and discussed in the paragraph below Eq. (|83[) . each such 
thick line segment can be thought of as a pair of coun- 
tcrdirected thin lines, corresponding to the propagation 
of a particle-hole pair. The constituent particle and hole 
lines carry Keldysh and frequency indices that traverse 
P\ without mixing. In Fig. [51 we picture the same basic 
diffuson propagator shown in Fig. [Jja), but with a twist 
of the right end relative to the left. The twist is rep- 
resented by the "oo" symbol. We will use this twisted 
representation of Pa whenever convenient to simplify the 
2D diagrammatic representation of the one loop correc- 
tions. 




FIG. 9: Basic diffuson propagator P A [Fig.[7£a) and Eq. (|82a|l 
with a twist. 




(b) 



(c) 







FIG. 10: Category 2)1: Diagrams renormalizing Xa- 



2. Renormalization of Xa 

Diagrams 231(a)-231(c) appearing in Fig. [TO] renor- 
malize the disorder parameter Xa- All three diagrams 
shown in this figure pair together two copies of the stiff- 
ness vertex V(a), using pairs of the disorder-only fast 
mode propagators Pa and Pa- There are no undeter- 
mined loop frequencies in these three diagrams, thanks, 
e.g., to the propagator twists in 231(a); the associated 
amplitudes therefore involve pure momentum shell inte- 
grations, since all propagator frequency indices are slow. 
[See Eq. ([75]) and Fig.0] We will compute 231(a) explic- 
itly; the labels 1 — * {tJi,a} and 3' — + {w 3 ,a'} in Fig. [TO] 
represent external, slow frequencies (u>i, W3) and Keldysh 
indices (a, a'). 

Using the Feynman rules, we find 



231(a) 



1 



d d q 



(89) 

where i and j denote vector components, with the pure 
momentum shell integration 



d 2 k 



(2tt)' 



(2k + q) l (-2k-q)^ 



X 



A£ a (| Wl |,M,k)A£ a (M> 3 |,k + q) 



dx 



2h 2 2D 2 (2n) Jx x 



-S^X 2 



2- 



-2dL 



(90) 



up to irrelevant slow mode frequency (0^1,0^3)- and mo- 
mentum (q)-dependent terms. On the last line of 
Eq. (|90p . we have performed the change of integration 
variables x = Dk 2 , with D = 1/Xh [Eq. (J73J)], and we 
have used Eq. ([55]). Combining Eqs. and using 
Eq. (|87[) . and summing over slow mode frequency and 
Keldysh indices, we obtain 



231(a) 



/ A 2 2dl \ X A 
V 8ttX a ) 2X 2 



d d r 



Tr 



(91) 



231(b) and 231(c) may be similarly evaluated; in fact, 
these diagrams exactly cancel because of the directional 
dependence of the stiffness vertex V(a) upon the loop 
momenta — see Table [H Thus the complete A^ renormal- 
ization is given by Eq. (JHU). 
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3. Renormalization of A and h 



The diagrams in category 2)2, Fig. [TTJ renormalize the 
disorder parameter A, proportional to the inverse dimen- 
sionless DC conductance. D2(a);-£)2(a)i v , shown in the 
left-hand column of Fig. 111! possess no undetermined 
loop frequencies, and therefore involve pure momentum 
shell integrations, similar to that in Eq. (|90| . In fact, 
S2(a)i and 2D2(a)iii exactly cancel D2(a)ii and £>2(a) iv , 
due to the presence of twists in the latter diagrams and 
the momentum-dependence of the stiffness vertex V(a) 
[Table IJ. 

On the other hand, the diagrams in the right-hand 
column of Fig. [TT] involve simultaneous frequency and 
momentum loop integrations, and their sum indeed yields 
a non- vanishing renormalization of A. 

Diagrams 3D2(b)j and 1)2(b)\\ pictured in Fig. [TT1 give 
identical contributions, and each involve two copies of the 
stiffness vertex V(a), with one basic diffuson and one 
interacting sector propagator P$ + Pc- We evaluate 



D2 ^ = ? ^^/^^^ a ^^ 

(92) 

where the frequency-momentum integral divides into two 
pieces 



case, up to irrelevant terms we have 
ij f du> a d 2 k 

2 25 



UJ a >0 



, . (2k) l (2k)Jr s (s Q - S2 )( S2 - Sa ) 
A a / (0, \w a | , k) A|' (K |,k)A°'' a '(|w a |,0, k) 

~2hT~ J l2^Ft k [Dk ~* " a] 

x[£>k 2 -i(l- 7s )r'a;a]" 1 }- 

(94) 



The factor (s Q — s 2 ) 2 appearing in the first line of Eq. (|94|) 
is inherited from Ps, which projects onto the propagator 
channel purely off-diagonal in sgn(w) space [see Eq. (|82c[) 
and the discussion below Eq. (|80b[) ] . Eq. yields a re- 
sult independent of the slow frequency w 2 > so we may take 
sgn(w 2 ) < without loss of generality. As a consequence, 
the only effect of the aforementioned "projection factor" 
in Eq. (|94p is the restriction of the loop frequency inte- 
gration to the half space uj a > 0, indicated on the third 
line of this equation. Following a change of variables, 
Eq. (94]) gives 



4tt 



2dl 



7s 



Ml ~ 7.) 



(95) 



T*,3 — Tl,3 

1 2 — 1 2S 



J 2C 



(93) 



corresponding to the Ps and Pc components of the in- 
teraction sector propagator, respectively. In the former 



The frequency-momentum shell integral J^{z;z ls ) in 
Eq. ([93)) is defined and evaluated in Appendix [E] 
[Eq. (|E6|)]. The CDW channel contribution jj£ to 
Eq. (|93p may be similarly computed, using the propaga- 
tor Pc, and working only to lowest order in r c [i.e. ignor- 
ing the logarithmic denominator, Eq. ([55)) . in Eq. (|82d|) ]. 
Using Eq. (|87p , and summing the identical contributions 
from 2)2(b)i and 3)2(b)n, one finds 



(a) 
(i) 

(ii) 

(iii) 

(iv) 









FIG. 11: Category 2)2: Diagrams renormalizing A. 



S2(b) i + S2(b) ii 



\2dl 



47T 



2 



x — / d rf rTr 
2A 



(96) 



If we set 7 c = in Eq. (|96jl , the remainder is recognized 
as the usual correction to the "dimensionless DC resis- 
tance" A in the presence of short-ranged density-density 
interparticle interactions ^ 62 ' 6 ^ 64 

The computation of the final two diagrams Q2(b)\\\ and 
£>2(b)iv in Fig.[TT]is more complicated, as each graph in- 
volves two loop frequency integrations. We will demon- 
strate the evaluation of diagrams involving two [lD6(a)ni] 
and three pD6(a) v ] frequency loops appearing in Fig. [TSI 



20 






FIG. 12: Category 333: Diagrams renormalizing h. 



below. Here, we simply give the result for the sum 



3?2(b)iii + 3?2(b)i. 



-X2dl 

47T 



2+— ^ln(l- 7s ) 



x — / d a rTr 
2 A 



(97) 



Obtaining this result requires the use of the integral for- 
mulae given by Eqs. (|E5[) and (|E13|) of Appendix [El We 
demonstrate shortly that the "anomalous" amplitude in 
Eq. (]9"T|) is precisely canceled by other diagrams. 

We turn now to the renormalization of h by the dia- 
grams in Fig. [T2l ©3(a) and 333(b) prove easy to eval- 
uate, involving only the associated symmetry-breaking 
vertex V(b), fused with the propagators Pa and P$ + Pc, 
respectively. 333(a) involves a pure momentum shell 
integration, while 233(b) requires the integral Ji(z,z'), 
Eq. JEj). The result is 



©3(a) + 333(b) 



-2dl 

~8^r 



{Xa + A [ln(l — 7 S ) + 7 C ]} 



x ih / cTrTr \w\&(Ql + Qs) 



Diagrams 333(c)-333(e) are as simple to evaluate, but the 
sum of the associated amplitudes gives zero. 

To complete the renormalization of A and h, we must 
compute the graphs pictured in Fig. [131 These diagrams 
represent frequency-momentum integrations quadrati- 
cally divergent in momentum; evaluation of the required 
integrals produces terms zeroth and first order in the 
energy cutoff A; the latter are presumably canceled by 
the measure, although we have not checked this in de- 
tail. 334(a) and 334(b) in Fig. rT3] involve the crosspair- 
ing of the interaction vertices V(c) and V(d), with both 
non-interacting Pa and interacting Ps + Pc propagator 
components. These diagrams are quite lengthy to evalu- 
ate, because they require two [334(a)] and three [334(b)] 
frequency loop integrations, and necessitate Taylor ex- 
pansions of the fast mode propagators in powers of the 



(a) 




(b) 





FIG. 13: Category 234: Diagrams renormalizing A and h. 



external frequencies and momenta. Other diagrams with 
multifrequency loop integrals will be tackled in detail in 
Sec. IIV A 5[ below. Here we merely quote the result 

334(a) + 334(b) = - 332(b) iM - 332(b) iv 
'\2dl 



8tt 

x / d d rTi 



[7 S + ln(l - 7s )] ih 



(Qs + Q s ) 



(99) 



In order to obtain Eq. ([99]) . one employs the frequency- 
momentum shell integral formulae given by Eqs. (|E10|) . 
(lElip . (|E14[) . and (|E15[1 . As promised, the "anomalous" 
correction to A obtained in Eq. (|97|) is completely can- 
celed by the amplitude 334(a) + 334(b). In addition, we 
pick up a crucial renormalization of h from Eq. 



4- Renormalization of the density of states v(uj) 

Before we treat the interparticle interaction parame- 
ters Tu and Ty, we pause to consider the local scaling 
operator 



dJ_ 
2tt 



-jf a,a 

•'UJ.Uj' 



(r) 



(100) 



With the aid of Eq. ([52"]) . it can be seen that the expec- 
tation value of v a {oj, r) in Eq. (|100[) represents a measure 
of the disorder-averaged, coarse-grained, single particle 
density of states (DOS) v{u>) in the diffusive Fermi liq- 
uid, with uj measured relative to the Fermi energy. [The 
integral over the auxiliary frequency oj' in Eq. (|100p is 
necessary to eliminate an energy-conserving delta func- 
tion in this expectation.] The scaling behavior of v(to) 
as a function of energy scale to or system size L may 
be determined through the renormalization of v a (uj,r). 
Applying the Wilsonian background field decomposition 
[Eqs. ([Ml) and ([71])] and the RG program [Eqs. (H]), 
([82)1 and Tabled] of the previous section to Eq. ([100]) , 
one encounters the same diagrams responsible for part 
of the renormalization of the frequency rescaling factor 
h, graphs 333(a) and 333(b), pictured in Fig. rTJ The 
category 334 "interaction-interaction" diagrams shown in 
Fig. 1131 which provide a further renormalization of h 
[Eq. (j99p ], do not appear in the computation of the scal- 
ing dimension of v a (cu,r). 
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We obtain the renormalization 



(a) 



(b) 



(^( W ,r)) A ~<^( W ,r)) x 

x (l + ^{AA + A[ln(l- 7s )+ 7c ]}) 

(101) 

similar to the diagrammatic amplitude expressed in 
Eq. JSgJ). In Eq. ([TOT]) , the symbols (. . .) A and (. . .) A 
denote the expectation value taken with respect to the 
generating functional Z in Eq. |76|) . before and after the 

elimination of the fast modes ^'^(k), respectively. In 
Sec. IIVB1 we will use Eq. (|101[) to derive a flow equation 
governing the scaling behavior of the DOS. 



Renormalization ofTu andVy 



(iii) 





d 






■ 


fa -a\ 


1 
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With the renormalization of the disorder-only sector 
parameters Xa, A, and h complete, we now turn to the 
(much more involved) renormalization of the interpar- 
ticle interaction parameters Tjj and Ty. It will prove 
convenient to introduce a set of six slow mode opera- 
tors {Ou,Ou,Oy,Ox,Ox,Oy}, defined by the expres- 
sions given m Table HH In this table, we have introduced 
the following compact notation for the position and fre- 



TABLE II: Slow mode operators that appear in the renor- 
malization of the interaction parameters IV an d IV (° r 
equivalently, T s and r c ). The operators in this table are 
summed over Keldysh species index (a), integrated over 
position space (r), and integrated over frequency indices 
{1, 2, 3, 4} — » {uji , u>2, <^3, ^4}; we have used the compact nota- 
tion introduced in Eq. (|102[1 for this integration. In the inter- 
acting sector of the pure slow mode FNLerM action, Eq. ((78} , 
the sum Ou + Ou couples to the same-sublattice interaction 
strength IV, while Ov couples to the intersublattice interac- 
tion constant IV. The operators Ox, Ox, and Oy do not 
correspond to terms occuring in the original FNLerM action 
[Eq. (|67[1 . after the saddle point shift in Eq. (|65[l ], but are gen- 
erated at intermediate steps in the renormalization process. 



FIG. 14: Category 2)5: Diagrams renormalizing Ft 
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quency space integral 
dcoi du>2 du>3 CL04 



d d v5 l+3 , 2+A = J d d v. (102) 
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Operators Ojj + Ou and Oy correspond to terms in the 
original FNLerM action, Eq. (|rI7|) [after the saddle point 
shift in Eq. (f65|) . with the replacement Q — > 1 + 5Q S , as 
in Eq. (|78|) ]. coupling to the same-sublattice and inter- 
sublattice interaction strengths Try and IV, respectively. 
The last three operators listed in Table UT1 are new; they 
do not correspond to terms appearing in Eq. ([rT7|) . The 
structures possessed by the operators Ox, Ox, and Oy 
closely parallel those of Ou, Ou, and Oy, with a cru- 
cial difference: the former lack the frequency-dependent 
vertex factors S1S3 — > sgn(uq) sgn(o> 3 ), etc., appearing in 
the latter. We will see that Ox, Ox, and Oy are gener- 
ated at intermediate steps in the RG process; in order for 
the theory to be renormalizable, these terms must com- 
pletely cancel upon summing all one-loop diagrammatic 
amplitudes. We will show that this indeed occurs^ 

We consider first diagrams that (nominally) renormal- 
ize the same-sublattice interaction strength IV, corre- 
sponding to operators Ou and Oy in Table llll We begin 
with the category 2)5 diagrams pictured in Fig. [T?] The 
elements of these diagrams include single copies of the 
V(g) and V(h) vertices, and one of each "flavor" of the 
propagator components Pa, Pa, and Ps + Pc- Diagrams 
235(a)i and 3D5(a)jj yield simple pure momentum shell in- 
tegrations, with the result 



S5(a)j +S5(a) n = 



(A - \ A )2dl 



4m 



TuO 



(103) 



£)5(a)jjj involves an undetermined loop frequency; we 
evaluate this diagram explicitly. Note that in provid- 
ing the frequency labels for the single fermion lines in 



22 



J)5(a)iii, Fig.Qj] we have neglected irrelevant small shifts 
of the loop frequency a — ► uj a by the slow mode fre- 
quencies {1,3} — > {cji,^}, necessary for strict energy 
conservation. We have 



S5(a) Mi = ^22-^+3,2+4^ 



d d r5Q^ 2 SQ 



'S3,4 J 3, 



(104) 

where the frequency-momentum integral divides into the 
pieces 



-f 3 = -f 3 S + ^3 C ■ 

We calculate /3 c; up to irrelevant terms, 
dw a <i 2 k 



(105) 



he = 



(2tt) 3 
x [A£ a (0, K|, 



(~s Q )r c (si + s a )(s 3 - s a ) 



(106) 



Pc yields the projection factor (si + s Q )(s3 — s Q ) in 
Eq. (|106|) : this factor gives a non-zero contribution only 
for s\ = —S3 = s a . We therefore obtain 



'3C 



r c (si - S3) 2 



h 2 



r c (si - s 3 ) : 



(2tt) 3 

_ \i£ a 2dl 
~ 32tt 



7c (si 



s 3 J 



(107) 



where we have used Eqs. (|85p and (|E3|) . Evaluating i 3 s 
[Eq. (|105[) ], in a similar fashion using Eq. (|E4[) . and com- 
bining this with Eqs. (|104[) and (|107|) . we find the result 



255(a) ni = 



A2 ^) [7c-Mi-7.)]r» (Ou - Ox) (108) 
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The amplitude in Eq. (|108[) renormalizes the same sub- 
lattice interaction operator Ojj «-> IV, but also generates 
its "evil twin," Ox, which does not appear in the original 
theory [Eqs. flBBJ) and (jrJTj) ]. 

The amplitudes for 255(b)i-255(b)ni mirror those found 
in Eqs. (fTU3| and !p8j) : 



255(b), + 255(b) M 



(A — A^)2rfl 
4^ 



(109) 



and 



DSCbiii, = ( ^) [ 7c - ln(l - 7 ,)]r p (Oy 



(110) 

We next direct our attention to the large array of 
diagrams shown in Fig. (fl5)) . denoted as category 256. 
Diagrams in subcategory 256(a) [336(b)] pair together 
two copies of the vertex V(c) [V(d)]. Moreover, as 
suggested by the figure, diagrams 256(a) m and 256(b) m , 
with m G {i, . . . , vi}, give structurally similar results that 
renormalize same-sublattice [Eq. (|62|) ] local operators in- 

volving SQs'i t 2$Qs3,4 (£>u or Ox) and 5Q f s "^Qsm (® u 
or Ox), respectively. We therefore evaluate the subcat- 
egories 236(a) and 256(b) simultaneously. 
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FIG. 15: Category 256: Diagrams renormalizing Ti 



256(a)i and 256(b)i each possess a single frequency loop, 
which may be evaluated using Eq. (|E3[) , giving the results 



2>6(a)i 

and 
2)6(b)i 



-X2dl 
8tt 

x [Tu{j s + j c )Ou + T v ( ls - lc )O x ] , (111) 



-X2dl 



8n 

x [Tu( ls + lc )Ou + r y ( 7s - 7c)Ox] • (112) 



256(a)ii and 256(a)iii give identical contributions. We 
evaluate 256(a)jj; as an example of a diagram with two 
frequency loop integrations over a — > uj a and j3 — > up. 
Using the Fcynman rules, we have 

256(a) Mi = lM!l!<5 1+3!2+4 J d d rSQ a s l 2 5Q a s l 4 h, 

(H3) 

where the frequency-momentum integral again separates 
into 



h — hs + he, 



(114) 



corresponding to the Ps and Pc components of the in- 
teraction sector propagator in 256(a),ji, Fig. [151 We first 
compute /4 s; up to irrelevant terms one finds 



/4S = 



du) a diop d 2 h 



r s (si - s Q )(s a+/3 - sp)s 2 a 



(2tt)* 

x (TuS3 - ^vs a )(TuSf3 - T v s a+ p) 

x [Ag°(Q > | Wa |,k)] 2 A|(|a; | l k). (115) 
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Although this expression entails two frequency loop inte- 
grations, the propagator kernels A© and A$ in Eq. qil5[) 
depend exclusively upon a single frequency variable, w a . 
The integral over up yields only a kinematical factor 
of \uj a \: the projector (s a +p — sp) in Eq. (|1 15[) gives 
a non-zero contribution as a function of up only over 
a frequency segment of length |w a |; over this range, 
— FySa+p) = Sf32T s . The factor (si — s a ) lim- 
its the range of the u> a integration to the half line (as 
usual); one obtains 



obtained are 



'4S 



-411 
irh 3 



{Tusis 3 + T v ) 



u) a >0 



d(jj a d 2 k 
(2tt) 3 



uj a [Dk 2 -i£ a Lu a y 



[D^ 2 -i{l- ls )^ a u a ]-\ 



We use Eq. (|E5jl to evaluate the integral, and we find the 
result 



US 



(T uSl s 3 +r v )J 2 (^;(l- 7s )e) 
(r l/Sl s 3 + ry)[ 7s +ln(l-7 s )], (117) 



-4TJ 

nh 3 
-X2dl 

32tt 



What we have found here is well-known (see e.g. Ref. |4S 
from other Finkel'stein NLcMs: in a one-loop calcula- 
tion, a given diagram may involve n = {0, 1,2,.. .} un- 
determined frequency loops; however, for n > 2, n— 1 of 
the associated frequency integrations typically give sim- 
ple kinematical factors, and the final diagrammatic am- 
plitude can almost always be expressed in terms of an 
integration over a single loop frequency and a single loop 
momentum, as in Eq. (|116|) [but see the evaluation of 
S9(b)i, Eq. dHD, and £13, Eq. (fT5I|) below]. 

We can try to evaluate /4 c in a similar fashion, but 
the associated amplitude goes like Ar 2 [l + 0(T C )]. In 
our one- loop results [Eqs. (|161[) and (|166[> . below], we 
keep only terms to second homogeneous order in A and 
T c in the flow equations for the interaction constants 7 S 
and 7 C ; therefore, we neglect the contribution of Zjc- 

Combining Eqs. (|113[) and (|116[) . and accounting for 
the identical contribution from 556(a)n, we have 



S>6(a)H + S>6(a)Hi 
and similarly 

S)6(b) M + seo^i 



X2dl 
Air 



[ 7s +ln(l- 7s )] 



x {TuOu + Y v Ox) , (118) 



X2dl\ , . 

—) K + Ml -%)] 
x (TuOu + T v Ox) ■ (119) 



Diagrams 2)6(a)j v and D6(b) lv also involve two fre- 
quency loops, and may be similarly evaluated; the results 



S36(a)i. 



and 



D6(b)i, 



-X2dl 
x (Ou- 

-X2dl 
(Ou- 



Ox), 



• ln(l- 7 ,)]-r c -i 



(120) 



r s [ 7s + in(i 



7,)]-r c f 



o 



x 



(121) 



D6(a) v and D6(a) v \ provide examples of the most com- 
plicated type of diagram that appears to one loop, incor- 
porating two copies of the interaction propagator compo- 
nent sum Ps + Pc , sandwiched between a pair of vertices 
[here V(c) and V(c)] that originate from the interacting 
sector [Eq. ([57]) ] of the FNLctM. These graphs involve 



(116) three loop frequencies a 



/3 — > up, and p 



u„ 



We evaluate J)6(a) v explicitly, as it requires less work 
than S6(a) v i. Using the Feynman rules, we write 

26(a) v = iM!l\ +3! 2 +4 J d d v8Q a s l 2 8Q a s l A h, (122) 

where now the loop integral ^5 breaks into jour pieces, 
/ 5S / 5 S S +/ 5 S C + / 5 C S + / 5 C C . (123) 

The two propagators occuring in the fast mode loop of 
£>6(a) v appear at the "top" and "bottom" of the associ- 
ated diagram in Fig. 1151 Each propagator may represent 
either Ps or Pc- In Eq. (|123p . the raised index describes 
the character of the "top" propagator component, while 
the lower index describes the character of the "bottom" 
component, i.e. i 5 c s is the amplitude for "top" compo- 
nent to be Pc and the "bottom" component to be Ps. 
Fortunately in this example, i 5 s c , I 5 C S , and J 5 C C all give 
corrections of order Ar 2 , which we ignore. 
We compute J 5 S S as follows: 



r s - 



du> a d/up du p d 2 k 



r s40l - S a ){s a+l3 - Sp) 



x (S3 - s a )(s a+p ~ s p )(T u s p ~ T v s a+p ) 

x {T uS p - T v s a+ p) [A^ a (0, |« Q |, k)A|(|w a |, k) 



(124) 



Similar to the two frequency loop case, Eq. (|115[) , the 
integrals over up and u p give only an overall kinematic 
factor of w 2 , thanks to the projectors {s a +p — sp) and 
(s a +p — s P ) in the integrand of Eq. (|124p . The factors 
(si — s a ) and (S3 — s a ) in this integrand allow for a non- 
zero amplitude only when s± — S3 = — s a . As a result, 
we obtain 



1 5S — 



4r 4 

^ x s 
7T 2 /l 4 



(S1 + S3) 2 

du n ePk 



u o >0 



(2tt) 3 



W 2 [L>k 2 -Z^ Q ]~ 



(125) 
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which may be evaluated using Eq. (|E8|) : 
4F 4 

(si + s 3 ) 2 



55 ~7T 2 ^ 

-Xi£ a 2dl 



1287T 



x r. s 



2-7. 
1-7, 



+ 21n(l-7.) 



(126) 



Combining Eqs. (|122| and (|126|1 . the resulting ampli- 
tude for 2)6(a) v is 



2)6(a) v 



-\2dl 



8tt 

x (Ou+Ox) 
Similarly, we find 

D6(b) v = 

x (OV + Ox) 



2-7.s 
1-7, 



2-7, 
1-7, 



21n(l- 7s 



(127) 



21n(l- 7s 



(128) 



Finally we give the results for the last two diagrams in 
Fig.QSl 



2)6(a) v 



X2dl 
8tt 



7s 



2-7, 
1-7, 



+ 21n(l- 7s ) 



-r, 



ln(l - 7,,) +7 S 



7? 



x (Ou - O x ) . 



(129) 



and 
2)6(b) vi = 



\2dl 
8tt 



I -— ^ ) +21n(l- 7s ) 



1-7, 



ln(l- 7 ,)+7, + y 



(130) 



x (Ou - Ox) • 



Obtaining these results requires the use of Eqs. (|E8[) and 
(ET0|. 

The diagrams in category 2)7, Fig.|T6l nominally renor- 
malize the intersublattice interaction strength Tv, corre- 
sponding to the operator Oy in Table [II] The computa- 
tion of these graphs closely parallels the results we have 
obtained above in Eqs. (p3]) . (fH^ - (jTT2l) . pTHjl - timj) . 
and (|127p — (|130|) . We therefore summarize only the re- 
sults: 



47T 



(131) 



2»7(a)ii 



»7(a)iii 



X2dl 
8tt 



A2d/ 
8tt 



[ 7c +ln(l- 7s )]r v (Oy + Oy) , (132) 

[ry( 7 , - lc)O v - Tuils + 7c)Oy] , 

(133) 



(a) 



(i) 




(ii) 




(Hi) 





(b) 
(i) 

(ii) 

(iii) 

(iv) 

(v) 













FIG. 16: Category D7: Diagrams renormalizing Tv- 



2>7(b)i + !D7(b)n 



-X2dl 

47T 



[ 7s +ln(l- 7 ,)] 



x (T v Ov - T v Oy) , (134) 



£7(b) m = (h^j |r s [ 7s + in(i - 7 ,)] + r c | 

x (O v + O y ) , 



(135) 



2)7(b)i, 



and finally 



A2dZ 
8tt 
x (O v 



r 

Oy) 



7s 



D7(b) v = 



-X2dl 
8tt 

+r, 

X (Oy + Oy 



r. 



1-7, 



+ 21n(l-7.) 

(136) 



+ 21n(l- 7s ) 



ln(l- 7s )+7 s + 



(137) 



Note that in Eqs. (|132p - (|137[) . we obtain both the renor- 
malization of the operator Oy *-* Ty , as well as the gen- 
eration of the unwanted term Oy, just as we had the 
generation of Ox and Ox in the discussion of Ou and 
Ou renormalization, above. 

Additional diagrams renormalizing the interparticle in- 
teraction strengths appear in Figs. [T7l l2"2l these will be 
discussed in the next subsection. Let us first pause to 
take stock of the results so far obtained. We will sum 
the amplitudes of all three diagram categories evaluated 
in this subsection, using the category symbols 2)5, 2)6, 
and 2)7 to denote the associated sums. We express our 
results in terms of the interaction parameters r s and r c , 



25 



defined via Eq. (|64|) , as well as the relative versions 7s 
and 7c, introduced in Eq. (|55j) . We drop terms higher 
than first order in T c (or 7c ). A large number of non- 
trivial cancelations occur, and the results for D5 + D6 
and Q7 prove quite simple. 

Summing categories £>5 and 2)6, Eqs. (|103[) . (|108|>— 
(fTT2l . (HHD-miD, and (Mlhllll, we find 



£5 + 1)6 



-2d/ 



8tt 



-2gM 
~8^T 



2 (Aa - A) (r s + r c )J (Ou + Ou) 

2\T sl: 



I -Is 



[Ox + Ox) 



(138) 



Summing the diagrams in category 2)7, Eqs. (|13ip - (|137|) . 
we obtain 



£7 = 



-2dl 



8tt 



-2dl 
~8^~ 



2Ar s7s 



1 



7s 



(139) 



In the partial results given by these equations, observe 
the almost complete cancelation of all "junk" terms in- 
volving the operators Ox, Ox, and Oy , not present in 
the original FNLctM action [Eqs. and (|57J)], In fact, 
the remaining terms proportional to Ox + Ox and Oy 
in Eqs. (| 1 38|) and (|139p . respectively, exactly cancel, up 
to terms fourth order in the slow mode fluctuations <5Q S 
and <5Qg- To see this, we use Eqs. (|69]l and (|70|) to write 



SQ s + SQl = -5Q s 5QI 



(140) 



Repeated applications of Eq. (|140[) prove the identifica- 
tion 



+ + (SQ s 5Ql) 



(141) 



Thus, to our working order in the slow mode fields, the 
FNLctM appears (so far) to be renormalizable. 66 We must 
complete the one-loop calculation to verify this. 



6. (Further) renormalization ofV 3 andT c 

The interaction sector renormalizations described by 
the remaining diagrams in Figs. 1171 1221 are most com- 
pactly stated in terms of corrections to the linear com- 
binations r s and r c [Eq. ([64])], which couple to the slow 
mode operators 



O s = Ou + Ou + 



v , 



and 



(a) 
(i) 



(ii) 



(Hi) 



V 



(iv) 




(b) 
(i) 

(ii) 
(iii) 



(iv) 



(v) 



(vi) 













FIG. 17: Category 2)8: Diagrams renormalizing T s = ru t Tv ■ 



respectively [Eq. (|T8[) ] . We choose now to quicken our 
pace, providing calculational details only when substan- 
tially different from those presented in the previous sub- 
section. 

Category £8 diagrams are shown in Fig. [ITJ Diagrams 
in subcategory £8(a) involve vertices V(j) and V(k), with 
propagator components Pa and Ps + Pc- The sum of these 
four graphs gives 



338(a) = 



-2dl 
x {T s O s + T c O c ) 



{2X A + 2X [ 7c + ln(l- 7s )]} 



(143) 



(a) 
(i) 

(ii) 

(iii) 



(iv) 









(b) 
(i) 



(ii) 



(iii) 



(iv) 




a (3 








O c = Ou + Ou - o v , 



(142) FIG. 18: Category 2)9: Diagrams renormalizing T s 
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The diagrams in subcategory 238(b) each pair together 
two copies of the vertex V(e), with various combinations 
of the propagator components Pa and Ps + Pc- Graphs 
338(b)i-5?8(b) V j each encompass two frequency loops, 
and their evaluation proceeds along the lines of 236(a);;; 
[Eqs. pT5|) - ([TT5|) ]. Diagrams 238(b)v and 238(b)vi involve 



three frequency loops a la Eq. (|124[) . but their amplitudes 
precisely cancel. The sum of all six graphs in 238(b) gives 



238(b) 



X2dl 
8tt 



2r c7 2 o s . 



(144) 



The eight diagrams in category 2)9, pictured in Fig. [151 
represent only two discernable amplitudes. These dia- 
grams also involve two copies of the vertex V(e), with 
the Xa propagator component Pa in combination with 
Pa or Ps + Pc- The four graphs in subcategory 339(a) 
each give identical contributions, with the sum 



»9(a) = (^) r,A 



(145) 



Likewise, each of the four diagrams in subcategory 
239(b) gives the same amplitude. Each of these graphs 
involves two undetermined frequency loops. We evalu- 
ate 239(b); explicitly, because the structure of the two 
frequency loop integration is quite different from that of 
Eq. (|115[) . studied in the previous subsection. Using the 
Feynman rules, we write 



239(b), 



m a ) 

2! 



x 2 



-'1+3,2+4 



d d r 



(siS( 
x(s 3 l 



^S3A 



s 2 SQ 



S 1,2) 



(146) 



In this equation, note that both of the two V(e) vertices 
contribute only the portions of the associated vertex am- 



(a) 
(') 

(ii) 




(b) 
(') 




(iv) 




(iv) 



FIG. 19: Category 2)10: Diagrams renormalizing 

p IV+rV 



plitudes [Table Q] proportional to T s ; all other terms van- 
ish. The integral Jg separates into parts 



'6S 



'6C, 



(147) 



corresponding to interacting propagator components Ps 
and Pc, respectively. From the figure, and given the 
structure of Ps, Eq. (|82c|) . it can be seen that /6S = 0. 
On the other hand, I e q may be written 

f dw a dwpd 2 k 2 2 2 

J (2^)4 Tc ( 2sa > ^l 3 ' k 

A£°(|a; a |, | WQ |,k)A£ a (| W/3 |, | W/J |,k)] 2 . (148) 

Distinct from the two frequency loop integration previ- 
ously evaluated, Eq. (|115|) . the propagator kernels A<p in 
Eq. (|148p depend upon either of the two loop frequencies 
Lu a or top. Eq. (|148p should be understood as an integra- 
tion over the full frequency-momentum shell, Eq. (|75[) . 
pictured in Fig. [5] in the space (|u>q.|> |w/3|, -Dk 2 ). It may 
be evaluated using Eq. (|E12|) : 



^6C 



2 6 r c 



J9 (2D 



~X 2 2dL 
ir(Trh) 2 



rT r . 



(149) 



Summing identical contributions from all four diagrams 
in 239(b), we obtain 



239(b) 



-X A 2dl 



2r c7 2 o s 



(150) 



The eight category 2310 diagrams depicted in Fig. 
pair V(e) with the other "triangular" interaction vertices, 
V(c) and V(d). The evaluation of these graphs proceeds 
as in Sec. IIV A 51 we give only the result for the entire 
category: 



2310 



(-^)r a {A[27c-21n(l-7.)] 



AXa1c)o s . 
(151) 



Categories 2)11 and 2312 appear in Figs. 1201 and [2T1 
respectively. The diagrams in these categories share the 
same structural elements as those in 2310, but there is no 
net contribution to the RG from either 2311 or 2312. The 
graphs in Fig. 1201 cancel pairwise: 



2311(a), + 2311(a) 



= 3>ll(a)n 
= 2311(b); 
= 2311(b);, 
= 0. 



^ 2311(a);, 

2311(b);; 

FSll(b)i, 



(152) 



The individual diagrams in category 2312 correspond to 
the generation of new operators, not present in the orig- 
inal FNLcrM action [Eqs. {66]) and {67])]. Fortunately, 
their sum gives zero [up to terms 0(6Q\), whose cance- 
lation we have not checked in detail] . 

Finally, we evaluate the single diagram in 2313, Fig. [221 
In this figure, two copies of the vertex V(f ) are connected 



27 



(a) 
(') 

(ii) 

(iii) 

(iv) 





(b) 
(') 










(iii) 



(iv) 







FIG. 20: Category 2)11: Diagrams renormalizing 



by a single, basic diffuson propagator P^. 3D 13 is unique 
among the graphs presented in this paper, in that it in- 
volves a pure frequency loop integration, since the mo- 
mentum carried by Pa between the slow mode fields at 
the vertices is necessarily slow. We will see that J) 13 
drives the CDW instability in the clean limit, A, — > 0. 
Such a contribution appears naturally in the frequency- 
momentum shell RG, but not in a pure momentum shell 
scheme^ i 3 ' 59 i 60 Using the Feynman rules, we write 



a\2 



£13 



(20 
2! 

x 



-x(s 3 5Q 



a, a 
S3 A 



S40Qs34 



(153) 



In writing the amplitude for 3D 13, Eq. (|153[) . we have ex- 
actly the opposite of the situation we had with 3D9(b)i, 
Eq. (fT46l) : both of the two V(f) vertices in Eq. (fl53")) con- 
tribute only the portions of the associated vertex ampli- 
tudes [Tablefl] proportional to F e ; all other terms vanish. 



(c) 



(b) 




(d) 





FIG. 21: Category 2)12: Additional non-zero diagrams, which 
cancel in pairs. Contributions from individual diagrams cor- 
respond to the generation of local operators not present in 
the original model. 
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Diagram renormalizing 



The pure frequency integral I7 may be evaluated as 
h=J ^(2 SQ ) 2 A£ Q (K|,| WQ |,q) 
■ dw a [Dq 2 - 2i£ a uj a ] ~ 

J A 



irh . 

4«r 



2nh 



2dl. 



(154) 



Combining Eqs. (|153|) and (] 1 54[) . we obtain the result 

5)13 = T clc 2dlO c . (155) 

This expression is second order in T c , but zeroth order in 
the disorder strength A. 



B. Dimensional analysis 

We obtain the flow equations for the FNLcrM parame- 
ters in the usual way, by re-exponentiating the diagram- 
matic corrections derived in the previous subsection and 
subtracting them from the pure slow mode action Ss, 
Eq. ([75]) . and by rescaling position r and time t via 



br, 

b z t, 



(156a) 
(156b) 



where few 1 + dl is the change in length scale, and z 
is the so far undetermined dynamic critical exponent. 
Equivalently, Eq. (|156[) implies that momentum k carries 
the engineering dimension 



[k] = 1, 



(157) 



while frequency to carries the (possibly scale dependent) 
dimension 



(158) 



both stated in inverse length units. Similar analysis of 
the FNLcrM action, Eqs. (55]) and (55]) [or §jo§ and (jBTj) ]. 
gives the engineering dimensions of the field matrix ele- 
ments 

= [Q a J(r)}=z, (159) 

and of the model coupling constants, 

[A] = [X A ] =2-d, (160a) 



[h] = [T s ] = [r c ] = d - z, 



(160b) 



28 



and 



[7s] = [ic] = o, 



(160c) 



where T s and r c ( 7s and 7c ) were denned in Eq. 
[Eq. (US])]. 

Combining the above dimensional analysis with the di- 
agrammatic results, Eqs. JSTJ, ([96 ]) — <[99 )i . (fl38|> . (fT39]) . 
pi5)) - (fTi5)) . (dnni), (HSU), and (HSH]), we obtain the one 
loop RG equations 



dX _ _ _ A^ 2A 2 r 
dj ~ 4^ 7c+ 4^ 

+ 0(A 2 7 2 ,A 3 ), 
dA^ A 2 Aa dA 



1 



■ln(l-7.) 



0(A 3 ), 



(161a) 
(161b) 



^ =(d - z)r s + i- (2A A - A) (r s + r c ) 



l 

47T 

0(Ar c7c , A^r^c, A 2 r s/c , AA^r 



2 (A- A A )(2- 7s )r s7c - A A r s7s 
/o AA^Tg/c), 



(161c) 



appears in the RG equations. The second order pure 
CDW self-interaction, (— 2r c7c ) in Eq. (|161d|) . is the sole 
exception to the rule; this term is tied to the ballistic 
Fermi liquid physics in the presence of sublattice sym- 
metry (and therefore nesting), and must be retained in 
the weakly disordered limit. We have checked by nu- 
merical integration that all of the results presented in 
Sec. El including the identification of a disorder-driven, 
interaction-mediated instability of the diffusive Fermi liq- 
uid phase in d = 2+e dimensions, are un-modified (within 
the perturbatively accessible coupling strength regime) 
by the inclusion of higher order terms in r c . 

Using Eq. (jlOip . one may also obtain a flow equation 
governing the scaling behavior of the disorder-averaged 
DOS v{u>). In the (Keldysh) sigma model formalism, 
Eqs. (fTUTJ)) and ^55]) imply that the DOS carries zero 
engineering dimension, so that Eq. (|101[) implies the flow 
equation 



^=-^[AA + Aln(l- 7s ) + A 7c ]. (163) 



^ =(d - z)T c + i- {2X A - A) (r s + r c ) 



1 

+ — 

47T 



2Ar c ln(l- 7s 



2r c7c 



o(Ar c7c , A^r c7c , A 2 r s/C , AA A r s/c ), (i6id) 



d\nh lr 

—y- ={d - z) + — [X A + A( 7c - 7s ) 

+ 0(A 7 2 ,A 2 ,AAa). 



(161e) 



In Eqs. (|161ap and (|161b[) . we have set d = 2 + e ex- 
plictly. The reader may be puzzled by the right-hand side 
of Eq. (|161bp , which seems to imply the wrong engineer- 
ing dimension (e) for Xa- However, the full engineering 
dimension (-e) arises through the sum of the first term in 
this equation with the third, which is itself proportional 
to Eq. (| 1 6 1 a|) . This third term in Eq. (|161b|) originates 
from the definition of A^: in the non- interacting sector of 
the Keldysh FNLcrM action, Eq. ([55]). the local "dimer- 
ization" disorder 67 , operator 



(Tr[Qt( r )VQ(r)]) 



(162) 



couples not to Xa, but to the ratio Xa/X 2 . Application 
of the chain rule to this ratio leads to the third term on 
the right-hand side of Eq. (|161bj) . 

The flow equations (|161ap - (|161ep are given to the 
lowest non-trivial orders in both the "dimensionless 
DC resistance" parameter A and the CDW interaction 
strength r c (or 7c ). Specifically, we retain terms in 
Eq. pMa]) - (fT6Tc|) . and (fTl)Te|) only to first order in T c , 
i.e. te rms Q (A 2 7c ) [Eq. |T6Ta|) ] and terms 0(XT c ,X A T c ) 
[Eqs. (|161cp and (|161ep ]. In principle, since we make no 
assumption about the smallness of Xa, we should retain 
terms 0(A^r c ), with n E {1,2,3...}, but only n = 1 



RESULTS AND DISCUSSION 



In this section we summarize and interpret the results 
of the RG calculation, set up in Sec. lIIII and performed in 
Sec. El for the FinkePstein NLctM (FNLctM) originally 
defined by Eqs. ([5? ]) -(|5g ]) in Sec. El The FNLcrM action 
given by Eqs. (|58]) and ([59]) is parameterized by five cou- 
pling constants: A, proportional to the dimensionless DC 
resistance of the system; A^, which gives a second mea- 
sure of the hopping disorder (associated with quenched 
orientational fluctuations of bond strength dimerization, 
as discussed in Sees. Ill A~3l and ln A 4p ; h, which tracks the 
relative scaling of length and time in the theory; Tjj and 
Ty, which characterize the strengths of generic short- 
ranged same-sublattice and intersublattice interparticle 
interactions, respectively, in the coarse-grained FNLcrM 
description. Alternatively, we have introduced the inter- 
action parameters T s and T c [the sum and difference of 
and Ty, see Eq. (|64"])]. which couple to the square of 
the smooth and of the sublattice staggered local charge 
densities, respectively. As discussed below Eq. (IM1) . a 
staggered interaction r c < is expected to promote 
the CDW instability in the clean limit [i.e. we identify 
r c ~ W c , with W c defined by Eq. © for the Hubbard- 
like model given by Eq. |TJ]. The RG calculation has 
been performed with the aid of an epsilon expansion in 
d = 2 + e spatial dimensions, with 0<eCl. 

The one-loop RG flow equations for the coupling 
strengths A, A^, r s , r c , and h were obtained at the end of 
the previous section, in Eqs. (|L6Lb[) (|161d[) . We restate 
them below in a slightly more streamlined form. We will 
need the "relative interaction parameters" 7s and 7c , de- 
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fined in Eq. (|85|) , repeated here for convenience: 
4 „ 2 



7s 

7c 



nh 
4 



irh 
2 

nh 



(T 



(164a) 
(164b) 



Next, we perform a trivial rescaling of the coupling con- 
stants appearing in Eqs. {55]) and (|161|) . 



A — > AttX, Xa — > 47tAa, 



(165) 



after which the one-loop RG flow equations for the cou- 
plings A, Xa, 7s, 7c, and ft. take the following form in 
d — 2 + e dimensions 



^Aa 
dZ 

dl 

dj c 
dl 

d\nh 
dl 



eX-X 2 lc 



-2X 2 



1 + ln(l - ls ) 



eX A + X 2 + 2- 



7s 

. A 4 <iA 



(166a) 
(166b) 



A dl 

Aa(1 - 7s) (7« + 2 7c - 27s 7c) 

-A(l - 7s) (7s + 7c - 2 7 s7c) , (166c) 

Aa (7c + 2 7s ) - A ( 7 s +7c) 

+A [2 7c ln(l - 7s ) + 7s7c ] - 2 7c 2 , (166d) 

(d-z) + A A + A( 7c ~7 s ), (166e) 



where I is the logarithm of the spatial length scale. 
The parameter z in Eq. (|166e|) is the (as yet unde- 
termined, possibly scale-dependent) "dynamic critical 
exponent."— We have also calculated the scale depen- 
dence of the disorder-averaged, single particle density 
of states (DOS) v{u>). Implementing the rescaling of 
Eq. (|165p in Eq. (|163|) . the one-loop flow equation for 
v is 



dims 
dl 



[X A + Aln(l - 7s ) + A 7c 



(167) 



Flow equations (|166ap - (|166e[) and (|167p are given to 
the lowest non-trivial order in the parameters A and 7c , 
but contain contributions from Xa and 7s to all orders. 
This is an advantage of the Finkel'stein NLcM formula- 
tion, which provides 3 - a loop expansion that is controlled 
perturbatively by the (small) dimensionless resistance A, 
but which does not require the coupling strengths Xa or 
7s to be small. We will discuss the physics that emerges 
from an analysis of our results [Eq. (|166p ], in turn for 
d = 2 and d > 2 dimensions. 



A. Structure of the one-loop flow equations 

Before turning to such an analysis, we pause to con- 
sider the structure of Eqs. ([T5oa|l - ([TB5e|l and (fW|) . To 



gain an understanding of the competing mechanisms 
driving the RG flow, we must attempt to isolate the ef- 
fects of the various terms. Fortunately, the most impor- 
tant structures in these equations cither occur generically 
in the perturbative description of interacting, diffusive 
Fermi liquids,— i & 62 i 63 i 64 or can be tied specifically to the 
special sublattice symmetry (SLS) [Eq. ((2J)] and its ef- 
fects upon the well-understood limits of (i) disorder with 
vanishing interparticle interactions ) 9 ' 12 ! 13 and (ii) inter- 
actions with vanishing disorder. 23 ' 24 Interestingly, as we 
will sec below, the competition between the disorder and 
interaction effects gives nevertheless rise to completely 
new physics. 



1. Non-interacting limit 

We consider first the non-interacting limit. In d = 2 
spatial dimensions, the FNLerM [Eqs. ([5T ]) -(|5g ]) ] with 
Tjj = Ty = (7s = 7c = 0) and non-zero disorder cou- 
plings A and A^ describes an unusual low-dimensional 
delocalized state^ ' 12 ' 13 analogous in many ways to the 
type of critical stat o 41 ' 69 ' 70 that occurs at a continuous, 
three-dimensional ("Anderson") metal insulator transi- 
tion (MIT) of non- interacting diffusive electrons, in the 
absence of SLS. At such an ordinary Anderson transition, 
the single particle wave functions at the Fermi energy are 
extended, and thus capable of transporting charge, heat, 
etc. across the system, but these states are also very far 
from the plane waves of a clean Fermi gas, and exhibit so- 
called "multifractal" scaling 69,70 (due to the presence of 
the disorder). In d = 2, the extended single particle wave 
functions reside only at the band center (energy u> = 0) 
for arbitrarily weak disorder [i.e. wave functions are ex- 
ponentially localized, with a localization length that di- 
verges upon approaching the band center (ui — > O)]^^ 
so that the critical state appears only at half filling [a 
necessary but not sufficient condition for SLS, Eq. ©]. 
The critical state at the band center (uj = 0) described 
above characterizes the 2D non-interacting random hop- 
ping model defined by the Hamiltonian in Eqs. ([T]) and 
©, with U = V = O^dl This delocalized state turns 
into a somewhat more conventional metallic one in d > 2, 
but it is still distinguished by SLS (and broken TRI). 
[See also the discussion in the Introduction, Sec. II Bl es- 
pecially Figs. and[U] 

These conclusions follow from the fact that a variety 
of aspects of the class AIII non-interacting sigma model, 
argued in Sec. Inland in Refs. I9l ll2lll3i to capture the low- 
energy physics of the sublattice symmetric random hop- 
ping model lacking TRI, can be solvecU 3 - exactly in d = 2 
using conformal field theory techniques. In particular, 
in the absence of interactions, it is possible 13 to obtain 
the exact renormalization group equations (to all orders 
in A and Xa) in 2D for the disorder-only sector sigma 
model parameters A and Xa, and for the average density 
of states v, in a particular RG scheme. These are of the 
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following form 

ir = °' (168a) 

— - h(X) + 2— -, 

= /i(A) (168b) 
^ = -Aa/ 2 (A) - / 3 (A), (168c) 

with fi(X), i € {1,2,3} real analytic (RG scheme depen- 
dent) functions of A. The second term on the first line of 
the right-hand side of Eq. (|168b[) , which vanishes here via 
Eq. p68ap . arises through the chain rule applied to the 
definition of Xa in Eq. (|58|). as discussed below Eq. (|161|) 
in Sec. IIVB1 Using the results of our perturbative anal- 
ysis [Eqs. (fT66a|) , (|166bj) . and p67]) ] and setting e = 0, 
we find agreement with Eq. (1168|) . obtaining the lowest 
order terms in the expansions 

/i(A) = A 2 + 0(A 3 ), (169a) 
/a (A) = 1+0(A), (169b) 
/a (A) = 0+O(A 2 ). (169c) 

[The term in square brackets on the second line of 
Eq. (|166ap vanishes when ^ s — ^ 0]. The lowest order 
(one-loop) results in Eqs. (|168[) and (|169p are universal 
and were originally obtained via perturbative methods 
by Gade and Wegner in Ref. 0. Eq. ([168 al) states that 
the "dimensionless DC resistance" A is exactly marginal 
in d = 2; equivalently, the dimensionless conductance 
g scales classically, since all quantum corrections to it 
vanish. The exact marginality of A demonstrates that 
Anderson localization effects are 'disabled' by SLS at the 
band center in the 2D random hopping model. We can 
extend these results for the one-loop RG equations of 
the NLerM to d = 2 + e dimensions, by simply adding the 
appropriate "engineering" dimension terms, i.e. — eA, to 
Eq. (fTSSaj) . and -eA^ to the second line of Eq. (|168bp . 72 
[See Eq. (|160a|l .] Moreover, this procedure yield o 73 ' 74 
the exact result for the RG equation for A, which [from 
Eq. (|168aj) ] reads in d = 2 + e dimensions: dX/dl = — eA. 
We conclude from this field theory treatment of the sys- 
tem that, in the absence of interactions and in dimen- 
sions d > 2, the delocalized phase, characterized by a 
finite conductivity, persists even for 'strong' disorder; in 
particular, even the strongly disordered system exhibits 
no transition into a localized phase (as already mentioned 
in the Introduction). 

The behavior of the system in d = 2 is special in that 
the "dimerization" disorder— parameter Xa is (logarith- 
mically) driven to strong coupling at a rate determined 
by /i(A) [Eq. (1168bj) ]. The running A^ feeds into the av- 
erage density of states as shown in Eq. (|168c|) . This 
equation in turn implies the divergence of the low-energy 
DOS upon approaching the band center, taken to reside 
at uj = 0. Specifically, one finds that in 2D, 9 j 12 ' 13 

v(w) ~ -exp(-c|lnw| Q ), (170) 



in the limit as ui — > 0, where c = c(A) is a scale in- 
dependent constant and the exponent a is 1/2.— Note 
that the DOS divergence in Eq. (|170p that occurs in the 
2D random hopping model has nothing to do with Fermi 
surface van Hove singularities, which may appear only in 
the clean limit,— e.g. in the case of pure nearest neigh- 
bor hopping on the square lattice, where such a singular- 
ity at half filling gives a weaker, logarithmic divergence. 
Returning to the random hopping model, in d > 2 the 
low energy DOS is finite, but parametrically enhanced 
at the band center. By contrast, the density of states in 
a non-interacting ordinary diffusive metal (lacking SLS, 
and being in one of the three Wigner-Dyson classes) is 
typically not renormalized by the disorder.— 

Summarizing, we have the following picture of the non- 
interacting random hopping model, which we have ar- 
gued to be described by the NLctM in Eqs. (|57l) - ([59)) . 
with Tu = Ty — 0: a delocalized phase exists (at the 
band center, i.e. for half filling) in d = 2 and (as always, 
and trivially) in d > 2; in neither case does a transition 
into an Anderson insulating state occur as the disorder 
strength is increased from zero. The single particle DOS 
diverges strongly upon approaching the band center in 
d — 2, while it is parametrically enhanced in d > 2. [See 
Figs. [3] and g]in Sec. [LBJ] 



2. Interparticle interactions: diffusive Fermi liquid and 
clean Hubbard-like model physics 

We turn now to the interpretation of various pieces of 
Eq. p66[) involving the interparticle interactions. We 
consider first the term in square brackets on the sec- 
ond line of Eq. (|166a|l . This term can be recognized 
as the usual perturbative correction to the inverse con- 
ductance in a diffusive Fermi liquid, in the presence of 
short-ranged interparticle interactions ^ 62 ' 63 and may be 
interpreted 6 ^ as coherent backscattering of carriers off of 
disorder-induced Friedel oscillations in the background 
electronic charge density. Background density fluctua- 
tions become a source of on-site disorder in the presence 
of electron-electron interactions, so that we may attribute 
this non-trivial correction to "dynamic SLS breaking." 

Next we note the non-trivial zero of the one-loop RG 
Eq. (|L66c|) at j s = 1. [A factor of (1— j s ) is expected in all 
orders.] This zero follows from the established represen- 
tation of the thermodynamic compressibility dnjd\i = k 
in terms of the Finkel'stein model parameters* 3 - 

k ~ h(l - 7.) = h - -T s , (171) 

valid in the diffusive Fermi liquid regime, where we have 
used Eq. (|164[) and where we set z — d [see Eq. (|166e[) ]. 
Eq. (|17ip shows that the incompressible limit k — > is 
attained by sending "f s — > 1 (for finite h). From the 
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definition Eq. (|164[) . we have 



dl 



1 d 
hJl 



1 dn 



•7.) 



-(1-7.) 
dln/i 



din h 
dl 



dl 



(172) 



In an interacting, disordered normal metal, k re- 
ceives no divergent corrections^ so that the condition 
(djs/dVj^s = 1) = is satisfied automatically. In the 
advent of sublattice symmetry, however, n does renor- 
malize, so that the first term on the right-hand side of 
Eq. (|172p is typically non-zero. This can be seen from 
the non-interacting limit, where n is equivalent to the 
single particle density of states v\ as shown in Eq. (|170p 
for the 2D case, v{u>) is strongly renormalized upon ap- 
proaching the band center (to — > 0). Regardless, it is 
plausible to expect that in the limit 7 S — 1 and k — 0, 
the incompressibility of the diffusive Fermi liquid is pre- 
served under the RG flow, so that dn/dl(-f s = 1) = 0, 
and therefore {d^s/dVj^s — 1) = 0, as we have found in 
Eq. UBBcJ). 

Finally, we note that the last term on the second line 
of Eq. (fl~66d)) drives the CDW instability, which is a rem- 
nant of the clean Hubbard-like model (recall that in our 
conventions 7 C < signals this instability). 



B. Results (i): No metallic phase in d — 2 

At last we analyze our results, considering first the 
2D case. As discussed above in Sec. IV A 1[ the non- 
interacting NLcrM describes in d — 2 an unusual crit- 
ical, delocalized state of the random hopping model in 
Eqs. {1} and ([9]) with U = V — 0. Other non-interacting 
models of disordered electrons are also known to possess 
delocalized states in d — 2, such as the spin-orbit ("sym- 
plectic") normal metal class, where e.g., a delocalized 
(metallic) phase occurs for sufficiently weak disorder^ 
The latter delocalized phase is however destabilized 3 
by the introduction of (in this case: long-ranged) 3D 
Coulomb interactions, yielding an insulator in the pres- 
ence of both disorder and interactions in d = 2. What is 
the analogous outcome in the presence of sublattice sym- 
metry, i.e. for the Hubbard-like models which are the 
focus of this work? Are the short-ranged interparticle in- 
teractions characterized by the (relative) parameters j s 
and 7 C relevant or irrelevant perturbations to the crit- 
ical, non-interacting system? The answer may be ob- 
tained by specializing Eqs. (|166a[) — (|166dj) to the vicinity 
of 7„ = 7 C = 0. The linearized flow equations are given 
by ipgap . (fl68b)) with pMaj) . and 



— ~ (3Xa 



2A) 7C/ , 



d\ v 
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-X A J\ 



(173a) 
(173b) 



where we have introduced the relative same-sublattice 
(intersublattice) interaction strength 7(7 = 7 S + 7c (7y = 



7 S - 7c) [c.f. Eqs. (64]) and (fT64)) ]. (Note that while this 
result is valid only to first order in A, no such restriction is 
placed on Xa-) Eq. (|173ap shows that the same-sublattice 
interaction strength 7(7 feeds back upon itself positively 
via the special sublattice class disorder coupling A ,4 > 0. 
Since the latter is always driven to strong coupling [see 
Eqs. (|168bj) and (fTMaj) ] in d = 2 for any A > 0, we see 
that same-sublattice interactions are rapidly enhanced 
as we renormalize the FNLcrM. At the very least, our 
result implies that the non-interacting description is un- 
suitable for describing the ground state of the 2D ver- 
sion of the full interacting, disordered Hubbard-like lat- 
tice model. This should be compared to an analogous 
result^ previously obtained for a TRI, interacting ran- 
dom hopping model on the honeycomb lattice. In both 
models, the non-interacting phase is initially destabilized 
by the growth of short-ranged same-sublattice interparti- 
cle interactions, as in Eq. (|173a[) . and this growth is me- 
diated by a special disorder coupling, here A^ , which oc- 
curs in general in the description of random lattice mod- 
els with an underlying SLS. We note that Eqs. (|173a|) 
and (|173b[) may alternatively be obtained from an anal- 
ysis performed in the (much simpler) non-interacting 2D 
(but disordered) theory (using e.g. replicas or supersym- 
metry) , since these equations describe only the scaling of 
the interaction operators near the non-interacting fixed 
point (i.e. are of linear order in the interaction strengths). 

Turning to an analysis of the full flow equations 
(|166ap - (|166ep , we observe that Eq. (|166b|) precludes the 
existence of any perturbatively accessible, non-trivial 
fixed points occuring at non-zero resistance A > 0, for 
small e > (d > 2). Numerically integrating Eq. (|166p 
for generic initial conditions in d = 2, we observe that 
either 7 C — > — 00, signaling CDW formation, or that 
A, Xa — > 00 and 7 C — > +00, indicating a flow toward 
simultaneously strong disorder and strong interactions. 
We demonstrate below that in d > 2 dimensions, these 
two flow directions away from the non-interacting state 
evolve into two distinguishable instabilities of the diffu- 
sive Fermi liquid, which exists as a stable phase through- 
out a region of finite volume in the four-dimensional 
(A, Xa, 7s, 7c) coupling constant space. By contrast, since 
we find no stable metallic region for e = 0, we expect that 
the 2D disordered and interacting Hubbard-like model in 
Eqs. (T]) and ([9]) is always an insulator at zero tempera- 
ture. 

We see that in 2D, sublattice symmetry is both the 
genesis of delocalized (critical) behavior in the absence 
of interactions, as well as the doom of such behavior in 
the presence of interactions. Moreover, in the Hubbard- 
like model studied in this work, SLS is also responsible 
for the Mott insulating charge density wave ground state 
in the clean limit. The physics that we have found is 
consistent with numerical studie s 39 ' 40 of the half-filled 
spin- 1/2 Hubbard model in d — 2, which have shown 
that TRI random hopping disorder preserves the charge 
compressibility gap of the clean Mott insulator, and that 
the disordered and interacting system shows no signs of 
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metallic behavior. 



C. Results (ii): Fermi liquid instability in d — 2 + e 

The situation in d — (2 + e) > 2 dimensions is more in- 
teresting. Upon increasing e from zero, a narrow, irregu- 
larly shaped sliver corresponding to a stable metallic, dif- 
fusive Fermi liquid state opens up in the four-dimensional 
(A, Xa, 7s, 7c) coupling constant space. [A metallic, diffu- 
sive Fermi liquid is a state characterized by the condition 
that both disorder parameters A and A^ tend to zero at 
large length scales, whereas the interaction parameters 
approach j s — » const, and 7 C — ► in the same limit; 
see Figs. [211 HU In this limit, the FNLcrM discussed in 
Section IIIII becomes a simple, weakly coupled Gaussian 
theory.—] The sliver encloses the line A = A^ = 7 C = 0, 
with — cxi < 7 S < 1, the entirety of which is perturba- 
tively accessible because the FNLcrM does not require the 
interaction strength ^ s to be small. A highly schematic 
3D "projected" phase diagram is depicted in Fig. [23] In 
this figure, the interaction constants (7 S , j c ) reside in the 
horizontal plane, while the vertical direction schemati- 
cally represents (both) disorder strengths, (A, Xa)- For 
example, we could take the vertical direction to measure 
the geometric average \JX Xa, while fixing the ratio X/Xa 
to a constant. The shaded sheath in Fig. [23] is a cartoon 
for the boundary of the stable metallic volume, which 
resides between it and the ballistic (A = A^ = 0) plane. 



"Anderson- 
Mott" 




/1 

7s 

FIG. 23: Schematic phase diagram in d = 2 + e dimensions, 
with < e <C 1. A and A a are both measures of the disorder, 
whereas 7 S and 7 C characterize the interaction strengths. The 
full phase diagram resides in the 4D (A,Aa,7s,7c) coupling 
constant space; this figure can be thought of as a 3D projec- 
tion of this space, where the vertical axis (perpendicular to 
the interaction axes) measures the total disorder strength for 
some fixed ratio of X/Xa- The stable metallic phase resides 
between the ballistic plane (A = Xa — 0) and the shaded 
sheath; the thick arrows indicate the two instabilities of the 
(metallic) diffusive Fermi liquid discussed in the text. 



Over the range of perturbatively small values of 7 C , the 
stable metallic, diffusive Fermi liquid phase resides in the 
region "f c > 0, terminating near 7 C = 0. For "f c > 0, the 
"height" of the stable metallic volume in the vertical "dis- 
order" direction is controlled by e, and is approximately 
independent of j c (as indicated in Fig. [23]). although the 
precise shape and size of the phase boundary does vary 
with the ratio X/Xa and with 7 S , making it difficult to 
characterize analytically. 

The flow equations (|166a[) - (|166d|) possess no nontriv- 
ial RG fixed points in dimensions d > 2, which would be 
accessible in an epsilon expansion about d — 2. Thus no 
continuous metal-insulator transition can be identified. 
However, the two instabilities described in the previous 
subsection for the d = 2 case persist for d > 2, and be- 
come clearly distinct roads out of the metallic state. The 
conventional CDW instability is found to always occur for 
initial "f c < and sufficiently weak disorder, i.e. when 
A, Xa "C e, and is represented by the flow j c — ► — oo. 
[Recall Eqs. (56]) and (fl~64"]) . and Section [LA]] This flow 
is accompanied by a decay in both disorder strengths 
A, A,4, indicative of the clean limit. In Fig. [23] the CDW 
instability is schematically indicated by the thick arrows 
that emanate from the terminating surface of the metal- 
lic phase near j c — 0, running off toward large negative 
7c- 

The primary result of this paper is the identification 
of a second route out of the diffusive Fermi liquid phase 
of the class AIII FNLcrM in d = (2 + e) > 2 dimensions, 
different from the Mott CDW instability, arising solely 
from the competition of disorder and interaction effects. 
As in d = 2 dimensions, this second route is charac- 
terized by a flow off to simultaneously strong disorder 
(A, Xa — > oo) and strong interactions (j c — > +oo), as 
indicated by the thick arrows emerging from the j c > 
portion of the phase boundary shown in Fig. 1231 we there- 
fore dub it an "Anderson-Mott" instability^ Even though 
there is no perturbatively accessible fixed point, this sec- 
ond instability is nonetheless perturbatively controlled 
in d = (2 + e) over a wide range of initial conditions 
when e << 1; in particular, it is accessible over the entire 
range < 7 S < IjZS Numerically integrating Eqs. (|166ap ~ 
(|166d[) for small e <C 1, we find that the "Anderson- 
Mott" instability^ can apparently always be reached by 
increasing only the "dimensionless resistance" A beyond 
some small threshold value At, while keeping the other 
three parameters A^, 7 S , 7 C fixed. The threshold value At 
is a function of e (and of the other coupling strengths) 
which vanishes continuously in the limit e — > 0. By con- 
trast, we find that it is difficult to access the "Anderson- 
Mott" route out of the Fermi liquid phase by varying the 
CDW interaction strength 7 C alone, despite the fact the 
ensuing instability is characterized by the rapid flow of 
7 C — * Too . We therefore interpret the boundary separat- 
ing the flow toward the stable metallic regime from that 
toward the regime of the "Anderson-Mott" instability as 
a disorder-driven, first order metal-insulator transition 
(MIT). We emphasize that a MIT does not exist in the 
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FIG. 24: RG flow trajectories in d = 2 + e dimensions, projected into the disorder (A,Aa) plane. In this and Fig. I25I 
e = 0.01. Each subfigure (a)-(d) depicts six different initial A values, distinguished by the labels 'A,'. . . ,'F', located at the 
initial conditions of the trajectories. The chosen values tune through the phase boundary enclosing the stable diffusive Fermi 
liquid state, demonstrating the disorder-driven instability discussed in the text. In each of these subfigures, the trajectories 
labeled 'E' and 'F' flow into the diffusive metal. As for the other coupling strengths, all trajectories in this figure share the 
initial conditions (ICs) Xa = — 7c = e = 0.01; trajectories in subfigures (a), (b), (c), and (d) share the IC 7 S = 0.1, 0.3, 0.5, 
and 0.7, respectively. The sequence (a)-(d) thus exhibits the 7 s -evolution of the phase boundary. 



non-interacting random hopping model, which possesses 
only a delocalized phase at half-filling for finite disorder 
in d > l , 8 ' 9 ' 12 i 13 while the clean spinless Hubbard model 
possesses only the Mott CDW instability. 

We demonstrate the physical picture described in the 
previous paragraph with a selection of RG flow trajec- 
tory plots. We have numerically integrated Eqs. (|166ap ~ 
(|166d[) for the case of d = 2 + e dimensions, with e = 0.01, 
for a variety of initial conditions lying within the pertur- 
batively accessible volume of the four-dimensional cou- 
pling constant space (A, \a,7s,7c)- Projected traces of 
these flow trajectories in the disorder (A, Xa) and interac- 
tion (7 S , 7 C ) planes are shown in Figs. l24l and l25l respec- 
tively. An individual trajectory is identified by a label A- 
F, drawn at the initial (final) condition of the correspond- 
ing projected trace in Fig. [2H (Fig. [23), an d by a subfig- 



ure label (a)-(d). Traces with matching labels in these 
two figures describe different projections of the same RG 
trajectory, which resides in the four-dimensional space 
(A, Xa, 7s, 7c)- In a given subfigure of Fig. [Ml six different 
initial values of the "dimensionless dc resistance" A have 
been chosen that straddle the phase boundary between 
the diffusive Fermi liquid state and the "Anderson-Mott" 
instability. In each of these subfigures, the trajectories 
labeled 'E' and 'F' flow into the metallic, diffusive Fermi 
liquid state ('diffusive metal'). As for the other coupling 
strengths, all trajectories in Figs. [M] and [25] share the 
initial conditions (ICs) A^ = — 7 C = e = 0.01; trajecto- 
ries in subfigures (a), (b), (c), and (d) share the IC 7 S 
= 0.1, 0.3, 0.5, and 0.7, respectively, as indicated by the 
black dots in Fig. [551 The subfigure sequence (a)-(d) 
thus exhibits the 7 s -evolution of the phase boundary. 
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FIG. 25: Same RG trajectories shown in Fig. 1241 here projected into the interaction (7s, 7c) plane. In each subfigure, the 
labels 'A,'. . . ,'F' denote [in contrast to Fig. [21] the endpoints of the trajectories bearing the same labels in the corresponding 
subfigure of Fig. 1241 Initial conditions (ICs) in the (7s, 7c) plane are shared by all trajectories in a given subfigure, indicated by 
the black dots. ICs in the disorder (A, Aa) plane may be read off from the plots in Fig. 1241 In each of the subfigures (a)-(d), 
the trajectories labeled 'E' and 'F' terminate in the diffusive metal. 



It is interesting to speculate upon the nature of the 
insulating state ultimately obtained upon breaching the 
boundary of the stable, metallic diffusive Fermi liquid 
along the "Anderson-Mott" route. Although we can- 
not rule out the possibility that both this and the CDW 
instability, while clearly distinct along the boundary of 
the diffusive Fermi liquid, eventually terminate into the 
same insulating phase, it is important to stress that 
this need not be the case. As shown in Figs. [M] and 
the "Anderson-Mott" instability is characterized by 
a rapid flow to large A a and large positive 7 C . Since 
Aa describes quenched orientational fluctuations in bond 
strength dimerization (as discussed in Sees. Ill A 31 and 
III A4[) . while j c is responsible for the Mott CDW insta- 
bility (7 C — > — oo), we might expect the insulator reached 
by the "Anderson-Mott" instability to possess local Mott 
insulating order, broken-up on longer length scales into 
randomly oriented domains pinned by the strong bond 
dimerization disorder. Note that a flow to large positive 
7 C is not consistent with local CDW formation [i.e. order 
at a nesting wavevector Kn, as in Eq. (01], but other 
types of Mott insulating order may occur in the limit 
of strong interparticle interactions (j c — > +oo). For ex- 



ample, sending the same sublattice interaction strength 
U — ► +00, for fixed V in the 2D square lattice version of 
the clean Hubbard- like model defined via Eq. (Q}, favors 
charge "stripes," i.e. order at K = (0, tt), rather than at 
Kn = (tt, tt). [Recall from Eqs. ([5f5)l and that this 
limit corresponds to 7 C — > +oo. In the same model, large 
negative values of 7 C would arise, e.g., from the different 
limit V — > +oo for fixed U ] 

Finally, let us put the "Anderson-Mott" instability, 
and the associated first order metal-insulator transition 
that we have found into some context. As we have 
stressed, the SLS possessed by the Hubbard-like model in 
Eqs. (P) and © leads to special properties in both, the 
non-interacting, but disordered limit (delocalized phase 
at band center in d > l)^ i 9 i 12 i 13 as well as in the inter- 
acting, but non-disordered (ballistic) limit (Mott insu- 
lating ground state at half-filling due to nesting for any 
V > U > 0). In the simultaneous presence of inter- 
actions and random hopping, we have a stable metallic 
phase, and have found an interaction-stabilized, disorder- 
driven MIT ind = 2 + e>2. These features mutate if 
we break SLS, e.g. either by allowing "diagonal" (on- 
site) disorder or by tuning the filling fraction away from 
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1/2. In the case of broken SLS (keeping broken TRI, 
for simplicity), we recover the usual unitary metal sym- 
metry class (for spinless electrons)^ described by an 
appropriate Finkel'stein NLcrM characterized by three 
coupling strengths, analogous to A, j a , and h appear- 
ing in Eqs. (l58|) and (|59|) . [Analogs of the "dimeriza- 
tion" disorder— parameter Xa, as well as the CDW in- 
teraction strength "f c do not appear in the unitary class 
Finkel'stein sigma model^ see the discussion in Sec. IIIII 
following Eq. (|86|) for details.] In the case of short-ranged 
interparticle interactions, the FNLuM predicts^ that the 
unitary metal class is an insulator in d = 2 regardless of 
the presence of interactions, while the MIT in d = 2 + e 
dimensions is actually in the same universality class as 
the non-interacting, continuous (pure Anderson) transi- 
tion, i.e. the interparticle interactions are irrelevant at 
the transition. (Recall that we consider the spinless case.) 



D. Conclusion 

As we have already summarized our main results in 
the Introduction (Sec. lIC|) . we conclude with a brief dis- 
cussion of applications and extensions of our work, and 
we comment on some very recent results^ for spin-1/2 
fermion systems subject to the effects of both disorder 
and interactions. 

As discussed in the last paragraph of Sec. IV C| the ad- 
dition of a SLS-breaking perturbation to the Hubbard- 
like model given by Eqs. |T]) and ([9]) is expected to 
push the system into (an interacting version of) the 
conventional, unitary metal class.— The (spinless) uni- 
tary metal with short-ranged interparticle interactions is 
known 3 - to exhibit a MIT in d = 2 + e dimensions, in the 
same universality class as the associated non-interacting 
(pure Anderson) transition in this symmetry class. It 
would be interesting to consider the crossover between 
the "Anderson-Mott" instability in the presence of SLS, 
identified in Sec. IV Cl and interpreted there as a disorder- 
driven, interaction-stabilized first order MIT, and the 
non-interacting, continuous unitary class Anderson tran- 
sition expected in the absence of SLS. 

Let us next revisit possible connections to spin-1/2 
fermion models. As discussed in Sec. Ill Bl and elabo- 
rated in Appendices [C] and [DJ the Finkel'stein NLcrM 
calculation presented in this paper also applies to (i) a 
spin-1/2 Hubbard lattice model with SLS, broken TRI, 
and strong spin-orbit coupling, and to (ii) a continuum 
model of fermionic quasiparticles, native to a supercon- 
ducting phase, possessing TRI and a remnant U(l) of 
spin SU(2) rotational symmetry (i.e. invariance under ro- 
tations about one axis in spin space), preserved in every 
realization of the disorder (polar p-wave superconduc- 
tor - see Appendix D). In the absence of interactions, 
both systems are realizations of the (quantum disorder) 
symmetry class AIII, within the classification scheme of 
Ref. 0. 

Consider first the interpretation in terms of the 



spin-1/2 random (sublattice) hopping model mentioned 
above. Disorder may be realized, e.g., via the applica- 
tion of a random orbital magnetic field to the otherwise 
clean model. (As discussed in Sec. Ill Bl and Appendix [Cl 
a random Zeeman field, on the other hand, realizes the 
completely different disorder class C, while the combina- 
tion of random orbital and Zeeman fields gives the ordi- 
nary unitary metal symmetry class A.) The calculation 
presented in this paper also applies to such a spin-1/2 
random hopping model, subject to generic short-ranged 
interparticle interactions, and to strong (homogenous or 
random) spin-orbit coupling. Spin-orbit coupling sup- 
presses the hydrodynamic spin diffusion channel in the 
Finkel'stein NLcrM low energy field theory, without mod- 
ifying the symmetry class (AIII) of the disordered system 
in the absence of interactions. [Without spin-orbit cou- 
pling, an additional interaction term between the spin- 
densities would appear, which does not exist in the spin- 
less case, and which is therefore not included in our treat- 
ment here.] 

On the other hand, the superconductor quasiparticle 
interpretation of the AIII NLcrM, mentioned above, pro- 
vides a very different realization of our results. In par- 
ticular, the superconductor interpretation does not re- 
quire the disorder to take any special form (i.e. pure po- 
tential scattering will work), but requires only TRI and 
a remnant U(l) of spin SU(2) rotational symmetry to 
be preserved in every disorder realization. As discussed 
in Appendix D, such a situation "naturally" arises in 
the description of a spin-triplet, p-wave superconductor 
(lacking intrinsic spin-orbit coupling) in its polar, TRI 
phase i 55 ' 56 This tantalizing, but so far speculative con- 
nection could prove a more readily attainable connection 
to experiment. 

As already mentioned, very recently Dell'Anna 3 ^ has 
independently studied several universality classes of 
Finkel'stein NLcrMs, including realizations of both the 
particle-hole symmetric class C and the sublattice sym- 
metric class AIII for spin-1/2 fermions. In the AIII case, 
Dell'Anna includes the spin diffusion channel (i.e. as- 
sumes no spin-orbit coupling). We have compared our 
one-loop RG results Eqs. (|166p and (|167j) to his, and we 
find agreement if we suppress the interactions associated 
with the spin degrees of freedom by hand. Dell'Anna has 
also computed flow equations for a spin-1/2 realization 
of the BDI class, appropriate to the half- filled spinful 
Hubbard model subject to real random hopping disorder 
[preserving TRI, SLS, and spin SU(2) rotational symme- 
try in every realization] . This case may potentially allow 
a detailed comparison with Monte Carlo simulations, be- 
cause the lattice Hubbard model with real random hop- 
ping manages to evade the infamous "sign problem," and 
numerical results in 2D are already available ! 39 ' At the 
time of this writing, however, Dell'Anna 3 ^ has not pro- 
vided a detailed analysis of his flow equations for this 
universality class (BDI). 
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APPENDIX A: SCHWINGER-KELDYSH 
SYMMETRY STRUCTURE 

In this appendix, we examine the symmetry content of 
the lattice fermion Keldysh action in Eqs. (fTT|) and (fT2j) . 
We establish a link between the non-interacting sector 
[Eq. (fTTj) ; see also (fT6"j) and (fT7|) ] and the transformation 
introduced in Eq. (|26|) . justifying the promotion of the 
latter to a spatially-varying "Goldstone fluctuation" in 
Eq. (USD. 

For a given realization of the complex random hop- 
ping amplitudes {6Uj}, the only "physical" invariance 
possessed by Eqs. (fTTl) and (fT2"j) is the discrete SLS, de- 
fined by Eq. (f^j). Nevertheless, the non-interacting sec- 
tor of the Keldysh action, Eq. (fTTj) . is actually invariant 
under a large group of continuous transformations, to 
be defined below. An analysis of the symmetry content 
of Eq. (jlip clarifies our derivation of the Keldysh sigma 
model provided in Sec. [TTJ in which the diffusive fluc- 
tuations of heat and charge in the non-interacting ran- 
dom hopping model emerge as Goldstone modes of this 
("spontaneously broken") continuous symmetry. The 
full, interacting Finkel'stein NLcrM, defined by Eqs. Q57[l- 
(|59p . further incorporates the interparticle interactions of 
Eq. (jT3J) , which formally obliterate this symmetry. 

Adopting the matrix notations introduced in the para- 
graph preceding Eq. (|14[) . we recast the non- interacting 
Keldysh action in Eqs. (fTS"j) and (fT7j) : 

Si ^Sf + Sf+S?, (Al) 

where 

Si 1 = ( CAi h,j c Bj + c B j t* j c A i j (A2a) 

(ij) 

S i = ir l^2 CAi^zizCAi +ir]^ CBj^isCBj, (A2b) 

i£A j£B 

and 

Si = 2jcAiS 3 |cj|cAi + 2J c B] t, 3 \uj\c Bj , (A2c) 

i€A j£B 

with t itj = t + St ip as in Eq. (fTTj) . 



Consider a general, spatially uniform unitary transfor- 
mation in |w| ® £ ® £ (frequencyCgiKeldysh) space: 

c Ai * Uj^CAi, CAi * C-AiV A , (A3) 

C B] U B C B j, C B j -> C Bj V B , 

where 

U\U A = V1V A = U B U B = V B V B = i H<9E{9e . (A4) 

We now examine the symmetry content of the non- 
interacting Keldysh action Si, Eq. (jAl[) . assembled in 
stages from the three pieces defined in Eqs. (|A2a[) - (|A2c[) . 

Invariance of the "Hamiltonian" piece SJ*, Eq. (|A2ap . 
under the transformation in Eq. (|A3|) requires that 

V A = Ul V B = U\. (A5) 

If we imagine discretizing and truncating the continuum 
of allowed absolute energies to a set of n discrete values, 
\uj\ — > \uj\ a , with a € {l,...,n}, then we can say that 
the conditions in Eqs. (|A4j) and (|A5p imply the action 
Si is invariant under transformations belonging to the 
space U(4n)xU(4n). 

Next, we consider the "Hamiltonian" plus "pole pre- 
scription" pieces, Eqs. (|A2aP and (|A2bp . respectively. In- 
variance of Si + Si requires the restriction given by 
Eq. ([A5[> . as well as 

t 3 i 3 U A t 3 i 3 = U B (A6) 

Taken together, Eqs. (|A5|) and (|A6|) imply that the 
U(4n)xU(4n) symmetry of 5/ 1 is broken down to the 
diagonal subgroup U(4n) by the pole prescription term, 
S?. 

Finally, the full non-interacting action is assembled by 
adding the "energy" piece, defined by Eq. (|A2cp . Un- 
der the transformation given by Eq. (|A3[) . invariance of 
Si = Si + Si' + S 1 " requires the conditions provided in 
Eqs. (fA5j) . (TATjj) . as well as 

[kU A ]=0, (A7) 

and 

[\u;\,U A }=0, (A8) 



with [A, B] = AB - BA. The non-interacting Keldysh 
action admits symmetry transformations off-diagonal in 
£ [sgn(w)] space, but necessarily diagonal in both \u\ and 
£ (Keldysh) spaces, implying the (much smaller) symme- 
try group [U(2) x U(2)] n . 

The pattern of symmetry breaking that we have found 
is summarized in Table [TTTl The symmetry content of the 
Hamiltonian piece, S/ 1 , of the Keldysh action is broken 
down into successively smaller remnants via the addition 
of the pole prescription and energy terms. We interpret 
the "pole prescription" piece S±, Eq. (|A2b|) with 77 — ^ + , 
as an infinitesimal source, analogous to an external field 
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TABLE III: Symmetry content of the non-interacting Keldysh 
action, Eqs. (|A1|) and (|A2|) . The target space of the NLcrM 
is set by the quotient of the first and second entries, below. 



Imposed condition 


Implied invariance 


1. Invc. of Si 1 : 


U(4n)xU(4n) 


2. Invc. of SJ* + S?: 


U(4n) 


3. Invc. of Si = + S? + Sf: 


[U(2) x U(2)] n 



in an O(N) ferromagnet; in the diffusive metallic phase, 
this source selects a particular ground state from the 
U(4n) xU(4n) saddle point manifold of the system. [The 
magnitude of the saddle point configuration is given by 
the elastic scattering rate 1/r due to the disorder, eval- 
uated in the self-consistent Born approximation. 41 See, 
e.g., Eqs. (HU and (|32")).] The (non-interacting) NLcrM 
retains only the quantum Goldstone fluctuations about 
this symmetry-broken minimum. On the other hand, the 
energy piece defined in Eq. (|A2c|) plays a decreasingly ef- 
fective role in suppressing fluctuations in the low energy 
limit, |cj| -C 1/r; we therefore retain modes (weakly) off- 
diagonal in both the \uj\ and Keldysh (£) spaces. As a re- 
sult, we expect that the non-interacting Keldysh NLcrM, 
appropriate to the complex random hopping model de- 
fined in Eqs. dJ) and ([9}, with U = V = 0, will possess 
the target manifold U(4n) x U(4n)/U(4n) - U(4n). The 
NLcrM with target manifold the Lie group U(4n), for 
finite, integer n, is termed the unitary principal chiral 
model in the field theory literature. An analogous re- 
sult for the non-interacting, class AM random hopping 
model, formulated in terms of bosonic replicas, was orig- 
inally obtained by Gade and Wegner^ 

Let us make several remarks. First, consider tuning 
the fermion density away from half filling in the non- 
interacting random hopping model. This necessitates the 
addition of a finite chemical potential term to the Hamil- 
tonian piece of the Keldysh action, Eq. (|A2a[) . Going 
through the above analysis, one finds a smaller NLcrM 
target manifold U(4n)/U(2n) x U(2n), characteristic of 
the TRI-broken, unitary normal metal class, as expected. 
We also note that the symmetry operation given by 
Eqs. (|A3[) and (|A4[) is not the most general unitary trans- 
formation that one may implement upon the independent 
Grassmann fields c? A , c£ A , c°g, and c°" B appearing in the 
Keldysh action [Eqs. (|A1|) and (|A2[) ]. In general, one 
should consider transformations that mix fields simul- 
taneously in sublattice flavor (c-a cb), particle-hole 
(cj «-» Ci), sgn(w) (E), Keldysh (£), and \w\ spaces. This is 
most easily done by concatenating all field degrees of free- 
dom at a particular location in position space into a sin- 
gle, many component Majorana fermion, upon which the 
largest unitary transformation possible is executed, and 
subsequently restricted in order to determine the sym- 
metry content of the Keldysh actionii 4 ^ This procedure 
is always required in the presence of TRI, in order to ob- 
tain the Cooperon modes. For the sublattice-symmetric, 



broken TRI (class AM) model studied here, only trans- 
formations diagonal in sublattice and particle-hole space 
preserve the Hamiltonian piece of the Keldysh action, S 1 / 1 
in Eq. (|X2a]) , 



APPENDIX B: GRADIENT EXPANSION 

We derive here the gradient expansion expressed in 
Eqs. (|48|) , (|49|) . and ([54]) for the fermionic functional de- 
terminant in Eq. (|39|) . Using the explicit construction 
[Eqs. (gU) and (|I5}] for the continuum NLcrM matrix 
field Q(r) in terms of its homogeneous saddle point value 
[Eq. (ST}], we have 



Q(r)=W(r) Qsptft(r), 



(Bl) 



where Qsp and U (r) possess the sublattice space decom- 
positions 



E 3 & 
£ 3 6 



and 



U{v) 



U A {v) 







U b (t) 



(B2) 



(B3) 



respectively. Note that Eqs. (|B1[) and (|B2[) imply that 

Q 2 (r) = i(l/2r) 2 . (B4) 

Now consider the determinant in Eq. (139|) . Abbreviat- 
ing Gq 1 + p = G _1 , we write 



S. 



DET 



Tr fin [G-i+i&MQsptiA ] 

Tr (in [ctiW (u^^G^U + iQ S p) ) 

Tr (in [i - i(2T) 2 ^a 1 G- 1 UQsp] ) 

Tr (in [-zQ S p(2r) 2 ] ) . (B5) 



Expanding the logarithm on the third line of Eq. (|B5|) . 
and dropping the constant term appearing on the fourth, 
we obtain 



where 



Si 



DET 



'DET 



i n {2rf 



s, 



(1) 

DET 



(2) 
DET 



■Tr 



(G^+p^ai 



(B6) 



(B7) 



Using the explicit form of the inverse Green's function 
G^ 1 in Eq. (|40| . as well as the sublattice decomposition 



of Q from Eq. 



the first term on the right-hand side 
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of Eq. (|B6| is 



4et ~«*r J ( (SaH + M7E3&) [Q f (r) + Q(r)] ) 

+ i2r J d d rTr ( M (r) Q(r) +p B (r) Qt( r ) j , 

(B8) 

where we have completed a partial trace over sublattice 
flavor (a) space. The first line of Eq. (|B8[) gives the term 
on the first line of the expansion in Eq. (|48|) . Writing 
out the trace on the second line of Eq. (IB8|) as a sum 
over the Keldysh species index {a} and an integral over 
time {t}, and using the locality of the auxiliary fields 
Pa/b P°a/b^) m tnese indices [Eq. |35|) ]. one obtains 
Eq. (HI. 

(2) 

The evaluation of the second term Sq ET in the expan- 
sion given by Eq. (|B6|) involves the terms in Eq. (|40|) that 
are functions of the momentum operator k, and are thus 
non-local in position space. Consider the components of 
Eq. (|40|) arising from a non-zero mean bond dimerization 
(<5^dim 7^ 0); Eq. ([25]) implies the microscopic definition 



r(k) 



i Stdim sin(k • n) 



(B9) 



(n is a unit vector specifying the orientation of the aver- 
age dimerization), so that 



i(j> I {-'k) = -i0j(k), 
(i^)*(k)=i^(-k). 



(BlOa) 
(BlOb) 



We define </>dim(r) as the following (coarse-grained, long 
wavelength) position space Fourier transform: 



^dim(r) = 



d d k 



e Jk - r *0/(k), 



(Bll) 



where the momentum integration should be taken up to 
some spherical cutoff of order the inverse mean free path. 
Eqs. (|B10p and (|B11[) imply that </>di m (r) obeys the con- 
ditions 



frdim(-r) = -^dim(r), 

#lim( r ) = 0dim(r). 



(B12a) 
(B12b) 



Thus </>dim(r) is a real, antisymmetric function. More- 
over, we expect 0dim( r ) to possess a "p-wave" (/ = 1 
angular momentum) structure as a function of r/r, given 
its momentum space representation, Eq. (|B9[) , which also 
implies that 



V0dim(r)| r =O 



(2n) d 



k^j(k) 



(B13) 



where c is some constant. 

At second order in the expansion given by Eq. (|B6 
Eq. JB7} gives 



Si 



(2) 
DET 



"(2r) 4 



Tr ([e(k) + ^(k) + &3*0/(k)] Q) 



(Tr[e(k)Q £ (k)Q 

e(k)QM<Mk)fi] ) ■ (B14) 



2Tr 



In Eq. (|B14|) . we have freely discarded terms of no rel- 
evance. In particular, we ignore ^(k), since this func- 
tion, although arising through the dimerization pertur- 
bation [Eqs. ([24]) and (|25[)]. is symmetric under reflec- 
tion, and does not generate anything not already present 
in the FNLcrM action (up to "trivial" anisotropy that 
may be removed via rescaling). Writing the first term 
in Eq. (|B14|) in position space {r}, expanding the band- 
structure e(k) ~ tTpk — + — ztJpV, with Up the average 
Fermi velocity (at half filling), and using Eq. (|44[) gives 
the usual "stiffness" term appearing on the second line 
of Eq. 



The second term in Eq. (|B14j) may be evaluated as 
follows: 
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2Tr e(k)g<7 3 i<Mk)Qj 

~t /" d d rTr ((-V 2 5 ri r')Q(r')^3^dim(r'-r)Q(r)) 



y d d rTr ((VV)-{2(r')^3 N dim (r'-r)] Q(r) + Q(r') <r 3 dim (r' - r) [vQ(r)]}) 
1 d d rTr [(-5 r y)[Q(r')a 3 [V 2 dim (r' - r)] Q(r) + Q(r') «r 3 [V0 dim (r' - r)] ■ [vQ(r 



+ Q(r') «r 3 [V0 dim (r' - r)] • [VQ(r)J + Q(r') a 3 dim (r' - r) 
2t y^rTr ( Q(r) «r 3 VQ(r) ) •VV dim (r') 
^ y d d rTr (Qt(r)VQ(r)) n, 



V 2 Q(r)]}) 



lr'=0 



CtStd 



(B15a) 
(B15b) 
(B15c) 

(B15d) 
(B15e) 
(B15f) 



where we have used an expansion for the bandstructure, 
e(k) oc tk 2 + 0(k 4 ), valid for small k (associated with 
the bottom of the band), t is the homogeneous hopping 
strength [Eq. |TJ)] , and where we have employed the com- 
pact notation £W(r — r') = 5 rr '. In going between sube- 
quations (|B15d[) and (|B15e[> . we have used Eqs. (IB4[) and 
(|B12a|> . while between Eqs. (|B15e|) and (|B15f|) we have 
employed the sublattice decomposition in Eq. (|44|) , com- 
pleted a partial trace over sublattice flavor space (a), and 
used Eq. (|B13|) . Combining Eq. (|B15f|l with the prefactor 
from (|B14[) . we obtain Eq. (|4"9"| (with a suitable redefini- 
tion of the constant c) . 



APPENDIX C: DISORDERED, BIPARTITE 
LATTICE MODELS FOR SPIN-1/2 ELECTRONS 

In this appendix, we identify the quantum disor- 
der symmetry classes appropriate to models of non- 
interacting spin-1/2 fermions, with nearest-neighbor hop- 
ping on a bipartite lattice at half-filling, subject to ran- 
dom orbital or random Zeeman magnetic fields. Our 
starting point is the general bipartite lattice Hamiltonian 



a,/3 



E c ^(^b)-/4+h.c. 

~ 2_^, C Ai \ L AA)i,i' c Ai' 



E4- 

3,3' 



(i 



BB>j,j> c Bj> 



(CI) 



where c A " and c B - are creation and annihilation opera- 
tors for fermions on the 'A' and sublattices, respec- 
tively, with a, (3 G {1, 1} spin-1/2 component indices. In 
Eq. (IC1|) . indices {«,?'} and {j,j'} respectively denote 
sites on the A and B sublattices, so that the sum on 



r 



runs over all (not just nearest neighbor) pairs of 
A and B sites, while the sums on and (J,j') run 

over all same sublattice pairs of sites. Particle hopping 
is facilitated by the spin-dependent intersublattice {t AB , 

t BB } complex connec- 



AB } and same sublattice {t AA . 



tivity matrices. The same sublattice matrices are con- 
strained by Hermiticity: 



t 

AA 



AA' 



t 



t 

BB 



BB- 



(C2) 



Let us adopt a compact matrix notation reminiscent of 
that implemented in Eq. P^|) . Sec. Ill A 11 We re- write 
Eq. jCU) as 



H = c r hc, 



(C3) 



where c (cr) is a column (row) vector, with indices 
in position, sublattice flavor, and spin-1/2 spaces, i.e. 
c — > c Ai ^ B j with all indices displayed. The single particle 

Hamiltonian h in Eq. 
decomposition 



C3j) has the sublattice flavor space 



h = 



^ AA ^AB 
ft f 
AB BB 



(C4) 



If we assume that the model in Eq. (|C1[) resides upon 
a e?-dimensional bipartite lattice of 2N sites, then h is a 
AN x AN Hermitian matrix. 

In the presence of disorder, the elements of h will 
be random variables. The structure of h can then 
be classified^^ using random matrix theory (RMT), 
according to the constraints imposed by any symme- 
tries preserved in every realization of the static disor- 
der. Consider the clean case of real, homogeneous, spin- 
independent nearest-neighbor hopping on a bipartite lat- 
tice at half- filling. We will refer to this situation as 
pure "sublattice" hopping. Pure sublattice hopping cor- 
responds to Eq. (|C1|) with 



t 



AA 



BB 



0, 



(C5) 
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and 

{i AB )t! =*%,->^> (C6) 

where t — t* , and 5u t j\ gives one for i and j nearest- 
neighbors, vanishing otherwise. We introduce two com- 
muting sets of Pauli matrices: the matrix a m acts in 
sublattice flavor space, while the matrix J n acts in the 
spin- 1/2 space, with m, n G {1,2,3}. We use the conven- 
tional basis for all Pauli matrices. In addition to symme- 
tries implied by the lattice space group, the single parti- 
cle Hamiltonian with pure sublattice hopping [Eqs. (|C4|) . 
(|C5|) . and (|C6|) ] possesses spin SU(2) rotational invari- 
ance 

Jm h Jm = h : (C7) 

for m — {1,2,3}, as well as three discrete symme- 
tries: time-reversal invariance (TRI), sublattice symme- 
try (SLS), and particle-hole symmetry (PH), defined as 
invariance under the following operations: 

c^-iJ 2 c, c+^c+iJa (TRI), (C8a) 

c-+a 3 {J) T , c^c T a 3 (SLS), (C8b) 

-U 2 a 3 (c t ) T , c f ^c T iJ 2 a 3 (PH). (C8c) 

In these equations, "T" denotes the matrix transpose op- 
eration. The TRI and SLS operations are antiunitary; 
the PH operation as defined in Eq. (|C8cp is a product of 
TRI and SLS, and is therefore unitary. The PH oper- 
ation so-defined is unconventional, in that it involves a 
spin flip. Eqs. (|C8a[) - (|C8c[) equivalently imply the fol- 
lowing conditions on h: 



J 2 h T J 2 = h 


(TRI), 


(C9a) 


-a 3 ha 3 = h 


(SLS), 


(C9b) 


J2&3 h T J 2 a 3 = h 


(PH). 


(C9c) 



Eq. (|C9bp demonstrates that SLS translates into a "cru- 
ral" condition on h: invariance under SLS forces h to be 
purely off-diagonal in sublattice flavor space (this condi- 
tion provides an equivalent definition of SLS in the ab- 
sence of interactions). 

For the case of disorder that breaks all of the "inter- 
nal" symmetries detailed in Eqs. (|C7|) and (IC9a|) (|C9c[) . 
the only constraint placed upon h is Hermiticity. 81 The 
most precise statement that we can make, independent 
of the details of a particular disorder realization, is that 
h belongs to a matrix representation^ of the Lie alge- 
bra u(4N) of the unitary group: we say that h G u(4N). 
This places the model in Eq. (|C1[) into the ordinary "uni- 
tary" (Wigner-Dyson) class AX Adding one or more of 
the symmetry constraints in Eqs. (|C7|) and (|C9a|) - (|C9c|) 
will allow us to refine this statement, potentially altering 
the random matrix class. 

Consider first the addition of a random orbital field to 
the pure sublattice hopping model defined by Eqs. (|C4|) . 



(|C5[) . and JC6J. The presence of a random orbital field 
modifies only Eq. (|C6[) . which now reads 

(i AB ):f^te^S {h3) S a >P, (CIO) 

with 8ij a bond-dependent real phase. Using Eqs. (|C5[) 
and (|C10|) in Eq. (|C4|) . we see that h preserves spin SU(2) 
rotational symmetry [Eq. (|C7|) ]. breaks TRI and PH 
[Eqs. (IC9al) and dCScl) ]. and preserves SLS [Eq. dC9bl) ]. 
Imposing Eq. (|C7|) alone implies that h G u(2N); subse- 
quently enforcing the "anti-SLS" constraint 

a 3 ha 3 = h (Cll) 

further reduces the space to which h belongs to u(N) x 
u(N). Since the conditions in Eqs. (|C9b[) and (|C11|) are 
complementary, enforcing instead spin SU(2) and SLS 
leads to the identification h G u(2N)/u(N) x u(N), asso- 
ciated with the "crural" random matrix class AIIL- Thus 
adding a random orbital magnetic field [Eq. (|C10[) ] to the 
spin-1/2 pure sublattice hopping model yields a system 
in the same class as the (non-interacting) spinless ran- 
dom hopping model [Eqs. |T]) and © with U = V = 0], 
discussed in Sees. II B l and IV AT] of this paper. 

We may also consider imposing only SLS, i.e. removing 
the spin SU(2) constraint given by Eq. (|C7[) . This is con- 
sistent with the addition of homogeneous spin-orbit cou- 
pling to the orbital magnetic field case analyzed above. 
Using Eq. dC9bl) . we find h G u(4N)/u(2N) x u(2N), so 
that h belongs to the same space of matrices (and there- 
fore the same random matrix class AIII), regardless of 
whether spin SU(2) rotational symmetry is preserved or 
destroyed (completely). As discussed in Sees. IIIBI and 
IVDl the Finkel'stein NLcrM formulated in Sec. [TTJ de- 
fined by Eqs. ([57| - ([59^) . also applies to a spin-1/2 Hub- 
bard model with sublattice hopping, subject to a random 
orbital magnetic field, and possessing spin-orbit coupling. 

Next we turn to the case of a random Zeeman 
field. Given the assumption of pure sublattice hopping 
[Eqs. (|C4|) . (|C5p . and (|C6[) ]. the introduction of a Zeeman 
field modifies Eq. (|C5|) as follows: 

R „ C12 

where J = {J m } — > {J^^}, with m G {1,2,3}, is a 
vector of spin-1/2 space Pauli matrices, and where b^i 
and hsj denote real three-component vectors with sub- 
lattice site-dependent orientations and magnitudes. Us- 
ing Eqs. (|C12|) and ([C6]) in ([Ci]) . we see that a ran- 
dom Zeeman field breaks spin SU(2) rotational symme- 
try [Eq. jCT}] (completely), TRI [Eq. dC9al) ]. and SLS 
[Eq. (|C9bp ]. while preserving PHH [Eq. (IC9cp ] in ev- 
ery static realization of disorder. Eq. (IC9c|) provides 
a symplectic condition on h; in isolation, PH therefore 
implies that h G sp(4N). The Lie algebra sp(4N) of 
the symplectic group is associated with the random ma- 
trix class CZ^l Class C is also realized by a spin-1/2 
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fermionic quasiparticle system native to a superconduc- 
tor with singlet pairing, broken TRI, and spin SU(2) 
rotational symmetry preserved in every realization of 

disorder 7,38,43,49.50,51.52 

In summary, given the assumption of a clean sys- 
tem of spin- 1/2 fcrmions with pure sublattice hopping 
[Eqs. ([Cl|, dHSl), and JC6|], the introduc tion of a ran- 
dom "orbital" magnetic field, as in Eq. (|C10|) . gives a 
non-interacting quantum disordered system with SLS, 
lacking TRI and PH, belonging to the random matrix 
class AIII; the introduction of a random Zeeman field, as 
in Eq. (|C12[) . gives a system with PH, lacking spin SU(2) 
rotational symmetry, TRI, and SLS, belonging to class 
C. The application of both random orbital and random 
Zeeman fields to a spin- 1/2 system with pure sublattice 
hopping breaks TRI, SLS, PH, and spin SU(2) rotational 
symmetry (completely), in which case the ordinary uni- 
tary (Wigner-Dyson) class A is recovered. 



APPENDIX D: POLAR SUPERCONDUCTOR 
INTERPRETATION 

Using the same type of random matrix theory (RMT) 
symmetry analysis 43 applied in Appendix [U] to classify 
systems of spin-1/2 lattice fermions, we demonstrate here 
a completely different "microscopic" interpretation of the 
FinkcPstein NLerM studied in this paper, which provides 
a view of our results alternative to the spinless Hub- 
bard model physics espoused in the main text. The 
microscopic view adopted in this appendix concerns the 
fermionic quasiparticles associated with a particular type 
of superconductor, subject to quenched disorder. (These 
quasiparticles are initially taken to be non-interacting, 
within a mean field theory treatment of pairing; subse- 
quently, quasiparticle interactions are added.) 

Fermionic quasiparticles in a superconductor do not 
carry well-defined quantities of physical electric charge.— 
In the presence of disorder, this means that electric 
charge is not a "slow," hydrodynamic variable; a low- 
energy non-linear sigma model description of such a dis- 
ordered non-interacting quasiparticle system can there- 
fore not describe any kind of conventional "metal- 
insulator" transition, as might be observed experimen- 
tally through electric transport measurements. Disor- 
der may yet dramatically influence other properties of 
the system, however, including thermal and spin trans- 
port coefficients [the latter only if at least a U(l) sub- 
group of spin SU(2) rotational symmetry is preserved], 
as well as the behavior of the low energy, single par- 
ticle (tunneling) density of states^ In particular, it is 
possible in principle to observe a "thermal metal" to 
"thermal insulator" transition ) 49 ' 50 ' 52 ' 85 indicative of An- 
derson localization of the single quasiparticle states due 
to the disorder. This transition is described by a ran- 
dom Bogoliubov-De-Gennes equation, and may possibly 
be accessible 5 ^ by tuning, e.g., the impurity concentra- 
tion in a superconducting sample. To build a more re- 



alistic model of such a system, one should also include 
the effects of interaction s 37 ' 38 ! 51 between the quasiparti- 
cles of the superconductor; near d = 2, both disorder and 
interaction effects may be reliably described within the 
Finkel'stein NLctM framework. 

In this appendix, we first show that a non-interacting, 
spin-1/2 quasiparticle system, native e.g. to a cer- 
tain spin-triplet p-wave superconductor and subject to 
quenched disorder satisfying certain additional symme- 
try constraints, falls into the same random matrix class 
(the "crural" class AIII in the scheme of Ref. [7]) as the 
non-interacting spinless random hopping model [Eqs. (TTJ 
and © with U = V = 0] discussed in Secs.lLBland lV A II 
To obtain a superconductor quasiparticle system in the 
"chiral" class AIII, the disorder must preserve, in every 
static realization, time- reversal invariance (TRI), as well 
as a remnant U(l) of the spin SU(2) rotational symmetry 
Such a situation is expected to occur "naturally" for pure 
potential scattering (i.e. non-magnetic impurities, and no 
spin-orbit scattering), in a spin-triplet p-wave supercon- 
ducting host residing in its TRI, "polar" phase i 55 ' 56 We 
emphasize that, in contrast to the Hubbard model, this 
quasiparticle system may be defined directly in the con- 
tinuum, without reference to a lattice or an additional 
sublattice symmetry. After demonstrating the equiva- 
lence of the non-interacting superconductor quasiparti- 
cle and spinless random hopping (normal) particle, or 
electron systems at the level of RMT, we will briefly dis- 
cuss the interpretation of the disorder and interaction pa- 
rameters of the corresponding FNLctM [Eqs. (|57p {59])], 
studied in this paper, in the superconductor quasiparticle 
context. 

In order to motivate the basic underlying idea, consider 
first the system of quasiparticles in a two-dimensional 
spinless (or: spin-polarized) p x -superconductor 
which is subject to static short-ranged disorder— (in the 
potential and pair field). Weak disorder, acting on the 
pair of Dirac Fermions at the two nodal points and pre- 
serving the TRI of the p x superconducting state^ 8 . is 
known 43 to place this system in the (orthogonal) chiral 7 
symmetry class BDL 89 ' 90 In the low-energy Dirac theory, 
the disorder occurs in two varieties: intra- and internode 
scattering. The intranode randomness 1 - 3 - takes the form 
of a random U(l) vector potential) 87 ' 91 In the p x super- 
conductor realization, quenched random U(l) vector po- 
tential fluctuations correspond to small random shifts of 
the positions of the nodal points, and thus, in particular, 
to small random fluctuations of the orientation of the 
Cooper pair wavefunction away from the cc-axis. Intern- 
ode scattering appears as a pair of random masses for the 
Dirac quasiparticles i 13 ' 87 In the disorder-averaged Dirac 
theory, it is necessary to specify two parameters <?a and g 
in order to quantify the strength of intra- and internode 
scattering, respectively. As shown in Ref. 13, <?a and g 
each play a role in the class BDI disordered Dirac theory 
described above that is analogous to the parameters Xa 
and A, respectively, in the NLctM formulation of the (dif- 
ferent) chiral disorder class AIII, as defined by Eqs. {57]), 
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(|58|) and (|60|) in Sec. HU We demonstrate below that 
this analogy reflects a concrete realization of class AIII 
in terms of quasiparticles of a certain p-wave supercon- 
ductor, with spin 1/2, rather than spin 0. The strength 
gA of the vector potential randomness in the spinless case 
is analogous to Xa in the sigma model description of the 
spin-1/2 quasiparticle system. Indeed, as we show be- 
low, Xa quantifies the strength of random fluctuations of 
the orientation of the Cooper pair wavefunction in the 
spin-1/2 p-wave superconductor discussed below. 

We turn now to the random matrix symmetry anal- 
ysis. Consider a clean p-wave, spin-triplet Bogoliubov- 
deGennes (pairing) Hamiltonian of the form 

H = f d d r C t(r) \-V 2 -fi] c(r) 



+ 



rd d r' 



b^r,!-') • A(r-r') +H.c. 



(Dl) 



where we define the spin-triplet pair field operator 

b(r,r') = c T (r) J 2 Jc(r'). (D2) 



In Eq. (jDip . c* and c denote two-component fermion 
creation and annihilation operators, i.e. c — > c a , where 
a € {Ti 1} is a spin-1/2 component index; fi is the chemi- 
cal potential. In the absence of disorder and quasiparticle 
interactions, the gap A(r) in Eq. (|D1|) is a static, spa- 
tially antisymmetric vector- valued function. In Eq. (|D2[) . 
the superscript "T" denotes the matrix transpose oper- 
ation, and we denote the vector of Pauli matrices act- 
ing in spin-1/2 space by J = {J m } — ► {Jm^}, with 
m G {1, 2, 3}. The pair field operator defined by Eq. (|D2|) 
carries electric charge 2e, and transforms like an 0(3) 
vector under SU(2)-spin rotations. In the "polar" phase 
of a spin-triplet, p-wave superfluid, the gap function takes 
the for m 55 ' 56 (ignoring dynamical fluctuations): 



A(r)=C0(r), 



(D3) 



where the fixed unit vector £ and the spatially antisym- 
metric scalar function </>(r) are both purely real. We take 
C = (0,0,1). The Hamiltonian in Eq. (fDTj) then pos- 
sesses time-reversal invariance (TRI) and a remnant U(l) 
of spin rotational symmetry (about the m — 3 axis), de- 
fined as invariance under the transformations 



-iJoc, 



cHJ 2 



[TRI], 
[spin U(l)]. 



(D4a) 
(D4b) 



The time-reversal operation [Eq. (|D4a[) ] is antiunitary. 

We implement the following compact notation: first, 
we construct the Majorana spinor 



X 



J 2 (ct) T 



(D5) 



which carries indices in particle- hole {p = 1,2}, spin-1/2 
{a =t, I}, and position space r, i.e. Xp=i( r ) = c a (r) and 



Xp=2( r ) = [^(c^) ] Q (r) with all indices displayed. Note 
that x transforms like c under spin space transformations 
[see e.g. Eq. (ID4b[) ]. Next, we introduce a second set 
standard of Pauli matrices {n m }, m £ {1,2,3}, acting 
in the particle-hole space of \- The definition Eq. (ID5P 
implies that 



X 



x»n 2 j 2 . 



(D6) 



Using Eqs. (ID5j) and (|D6j) . we re- write Eq. (|Dlj) as 



H = -x r h X , 



(D7) 



where we have introduced the single particle Hamiltonian 



h = 



J A 



A J 

-J2t T J 2 



(D8) 



with 



t -» ^(r, r') = [(-V 2 - f i)5 d (r - r')]S a ^ , 

(D9) 



A • J -» A(r - r') • J a ' P 



(a and [3 are spin-1/2 component indices.) 

Combining Eqs. (|D6[) and (|D7|) . and noting that the 
components of \ mutually anticommute up to irrelevant 
additive constants, we see that the single particle Hamil- 
tonian h may be taken without loss of generality to satisfy 
the "Majorana" condition 



n 2 J, h T J 2 n 2 



(D10) 



The transformations defined by Eqs. (|D4a[) and (|D4b|) 
may be re-expressed in the x language using Eq. (|D5jl : 
doing so shows that TRI and spin U(l) rotational sym- 
metry translate into the following conditions upon h [via 
Eq. jD3) and using Eq. (|DT0)) ]: 



Tii ft.ni =h 

JshJ?, =h 



[TRI], 
[spin U(l)]. 



(Dlla) 
(Dllb) 



Eq. (jDllap shows that time-reversal invariance appears 
as a "chiral" 7 condition [see Eq. (|C9bj on the Bogoliubov- 
deGennes single particle Hamiltonian h. 

As discussed in Appendix [Cl in the presence of disor- 
der, we may classify the non-interacting quasiparticle sys- 
tem introduced in Eq. (jDip by examining the algebraic 
structure of the Hamiltonian h [Eqs. (|D7|) and (|D8|> ] . aris- 
ing in the presence of an arbitrary, static realization of 
impurities. If we imagine discretizing position space into 
a lattice of N sites, then h is a AN x 4A^ Hermitian square 
matrix. In the absence of both TRI and spin U(l) rota- 
tional invariance, the only constraint upon this matrix is 
the "Majorana" condition, Eq. (|D10[) , which implies that 
h belongs to a matrix representation of the Lie algebra 
so(4N) [in the defining (vector) representation]. We write 



43 



this condition as h G so(4N). This space of matrices h is 
associated with the random matrix class D^ ' 43 i 85 

Let us turn to the situation of interest, a quasiparti- 
cle system, subject to (e.g. pure potential) disorder, that 
preserves the symmetry conditions native to the "polar" 
phase of a p-wave, spin-triplet superconducting hos t 55 ' 56 : 
TRI and spin U(l) rotational invariance. The "Majo- 
rana" [Eq. (|DlOj) ] and spin U(l) [Eq. (|Dllb[) ] conditions 
imply the following decomposition in spin-1/2 space: 



h = 



l 2N 






-n 2 h T 2N n 2 



(D12) 



where h 2N G u(2N). Further imposing the "anti-TRI" 
condition 



Hi /i IIi = h 



(D13) 



would identify h 2N G u(N) x u(N). Eqs. (|Dllap and 
(|D13[) are exact complements, so that instead, imposing 
the "Majorana" condition, spin U(l) rotational invari- 
ance, and TRI leads to h G u(2N)/u(N) x u(N). This 
space of matrices is associated with the "chiral" ran- 
dom matrix class AIIIil The same space of matrices also 
applies to the non-interacting, spinless random hopping 
model discussed in Sees. II B l and IV ATI of this paper 3^ 
[See also the discussion below Eq. (|C11[) in Appendix [Cl] 

At the level of random matrix theory, we have suc- 
ceeded in identifying (i) the system of non-interacting 
spin-1/2 quasiparticles (treated within the mean field 
theory of pairing) native to a superconducting host and 
subject to quenched disorder preserving TRI and a rem- 
nant U(l) of the SU(2) spin rotational symmetry, with 
(ii) the non-interacting limit of the spinless Hubbard- 
like model with random hopping, defined by Eqs. ([T]) 
and ([9]) with U = V = 0. We have shown how such 
a system might "naturally" arise in the context of a 
p-wave, spin-triplet superconductor in its TRI, "polar" 
phase ; 55 ' 56 subject to pure potential scattering due to 
non-magnetic impurities. 

We can extend the analogy between these two sys- 
tems to the detailed NLcrM formulation, incorporating 
residual quasiparticle interactions a la Finkel'steini^ In 
contrast to the 2D case of Dirac quasiparticles in the 
spinless p x superconductor mentioned in the paragraph 
above Eq. (|D1[) , we now restrict our discussion to the case 
d > 2, for which our RG calculation of Sec. IV CI predicts 
a metal-insulator transition due to the interplay of both 
disorder and interactions, as we will now explain. 

In Sec.|TTJ we derived the class AIII Finkerstein NLctM 
(FNLcrM), using the Hubbard- like model [Eqs. (TJ) and 
©] as our "microscopic" starting point. The action for 
the FNLcrM was defined in Sec. Ill A 41 by Eqs. (JSSJ) and 
(f59"| . In order to understand the relevance of Sec. IV CI 
in the superconductor context, we need to explain the 
re-interpretation of the FNLcrM parameters appropriate 
to the p-wave superconductor quasiparticle view, which 
we will now do^ The parameter A still plays the role 



of the "dimensionless resistance", proportional now to 
the U(l) spin-conductivity, associating positive and neg- 
ative spin-U(l) charges to particles with 'up' and 'down' 
m = 3-component of spin, respectively 38 ' 49 ' 50 The sec- 
ond disorder strength A^ was attributed in Sec. Ill A 41 
to quenched bond dimerization fluctuations in the ran- 
dom (sublattice) hopping model. As discussed in the 
paragraph above Eq. (|D1[) . in the superconductor \a 
measures the strength of quenched random orientational 
fluctuations of the p-wave Cooper pairing wavefunction 
induced by the disorder (the argument is outlined in foot- 
note |93|). 

Turning to the interaction sector of the FNLcrM action, 
Eq. (|59p . we now re- interpret the interaction strengths 
T s and r c [Eq. within the context of the quasipar- 
ticles of the superconductor. The parameter T s arises 
from the S z -S z component of spin-triplet* 9 - 4 - interactions 
(in the particle-hole channel) inherited from the normal 
Fermi liquid phase adjacent to the BCS superconductor, 
and modifies the effective spin diffusion constant in the 
presence of the interactions. 38 - The parameter r c , on the 
other hand, characterizes the residual interaction in the 
particle-particle Cooper channel, and may be interpreted 
for r c < as an attractive BCS interaction in a differ- 
ent, spin-singlet (e.g. s-wave) pairing channel^ Summa- 
rizing, the interaction sector of the Finkel'stein NLcrM 
defined by Eq. ([59")) can be re-expressed as follows: 



5/ =iJ2^ dtddr r 



")t a i a 



(r) 



h't 



t a i a I 



■ e a / dtd d v [r s 

a=l,2 ^ 



( r ) _ ^t.t 

S a z {t,v)] 2 
r c A£(i,r)A£(t,r)l 



(D14) 



where the fields 5f (t, r), A|(t, r), and Ag(i, r) may be 
written in terms of Q"' t a (r) and Q\l' a {r), and represent 
the electron bilinear operators 

S a z (t, r) ~ T [cf (t, r) dftt, r) - gj(t, r) c£(t, r)] , (D15) 
A§(«,r)~e°cJ(t,r)c?(t,r), (D16) 
AS(t,r)~e»e?(t,r)eJ(i,r). (D17) 

As written, A^(t, r) represents the simplest singlet chan- 
nel pairing amplitude, a G {1,2} is the Keldysh index, 
and £ a was defined by Eq. (fH 



APPENDIX E: FREQUENCY-MOMENTUM 
SHELL INTEGRALS 



The frequency-momentum loop integrations needed 
in the RG calculation described in Sec. IIVI are cata- 
loged in this appendix. All integrations listed below are 
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J 12 (z; z') =—d z J w (z; z) 
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tions enter (up to irrelevant perturbations), because easy- 
axis spin is conserved, in analogy to the situation with 
uniaxial spin symmetry discussed in Ref. [3^ 
The spin singlet nature of the residual pairing channel 
characterized by F c is encoded at the coarse level of the 
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